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TECHNICAL MEMORANDUM X- 64744 


COMPUTER USER'S MANUAL FOR A GENERALIZED CURVE 
FIT AND PLOTTING PROGRAM 

SECTION I. INTRODUCTION 


The use of graphically displayed data is often quite beneficial to engi- 
neers, scientists, and other decision makers. Such graphs can be in the 
form of reporting charts, presentation viewgraphs, and technical illustrations. 
A program has been developed for the IBM 1130 1 computing system (with a 
1627 plotter) that enables a non-computer oriented individual to easily gen- 
erate a plot of a set or sets of numerical data on 8. 5 x 11 in. standard 2 paper. 
This program is entitled the Generalized Curfit and Plotting (GECAP) program. 
This computer program has various output options in which it can display the 
data in different forms. These options include discrete, linear, continuous, 
and histogram types of output plots. The continuous option also performs a 
least-squares curve fit (Regression Analysis) on the input data. A statistical 
analysis of variance test determines the best-degree polynomial fit, and a 
printout of the polynomial coefficients is given. 

Graph titles, coordinate labels, and symbol identification are included 
on each plot. The user of the program needs only to supply two instruction 
cards and the necessary data to be plotted. 

The rest of this manual contains information about the use and opera- 
tion of GECAP. Five example plots and associated input data are illustrated. 
These examples demonstrate various combinations of the input options avail- 
able to the program. A brief summary of the program description is included 
along with Fortran listings of the coded subroutines. Modification of the pro- 
gram can be easily undertaken by experienced programmers since the logic 


1. This program can be readily converted over to the other computing sys- 
tems with minor modification. 

2. Any type of paper may be used; however, this size is most common in 
written reports and viewgraphs. 


operations within the program listing have been well illustrated. A compre- 
hensive list of various error messages is also contained within the program 
logic. 


SECTION II. CONTROL CARDS, INPUT DATA 
AND OPERATING INSTRUCTIONS 


A. Control Cards 

The following is a description of the control cards necessary for the 
program. 

• Control Card 1 may be either a cold start card or a job card 
A cold start card is mandatory each time a different disc is changed out of 
the computer. A job card has the following form: 

// JOB 

• Control Card 2 is an execution card and has the following form: 

// XEQ GECAP 1 

• Control Card 3 initiates the segmentation of the program in order 
to minimize the core requirements for GECAP. This card has the following 
form: 


*LOCALGECAP, LABEL, HIST, ERRO 

B. Input Data 

The following is a description of the data cards necessary for the 
program. 

• Case Card 

The format for this card is (A4, 12, IX, A2, 12, 3X, 5A6, 3X, 

12, 3X, 12) 
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Columns 1-4: Punch the letters CASE. 


Columns 5-6: Punch the number of the case. ( Punch 1 if the 
first case, 2 if second case, etc.) 

Columns 8-9: Punch the letters OF. 

Columns 10-11: Punch the total number of cases to be plotted 
per graph. 

Columns 15-44: These spaces are reserved for the name or 
label the user may wish to place on each individual data set. 3 
(Optional, leave blank if not needed) 

Columns 48-49: These columns are reserved for the exponent 
used to label the x axis. 4 (Optional, leave blank if not needed) 

Columns 53-54: These columns are reserved for the exponent 
used to label the y axis. 4 (Optional, leave blank if not needed) 

Columns 55-80: These columns should always be left blank. 

An example of this card is given as follows: 

Columns 

1-4 6, 8^9 11 15-44 48-49 53-54 

CASE 1 OF 4 TEST RUN 24D -1 -1 

• Program Instruction Card 

This card contains information about how the grid should be 
drawn and how the data should be plotted. The format for this card is 
(6F10.4, 4A1) 


3. This name will be printed below the x axis and will be preceded by a data 
symbol corresponding to the data mark symbols on the curve of interest. 

4. These exponents will cause the last labeled 'tic* mark on each axis to be 

followed by a 'x 10 eX ^ , where exp is the input value. (-9 s exp ^ 99). 
Example: 20 x 10 -1 . See Example Case 4 for use of this input variable. 
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Columns 1-10: XINC - This variable is the numerical incre- 
ment between ’TIC’ marks for the x axis. 


Columns 11-20: 
Columns 21-30: 
Columns 31-40: 
Columns 41-50: 
Columns 51-60: 


YINC - Same as XINC except for the y axis. 
XMAX - Upper limit for the x axis. 

YMAX - Upper limit for the y axis. 

XSTRT - Value assigned to the x origin. 
YSTRT - Value assigned to the y origin. 


NOTE: A restriction of the program requires that none of the 
six variables defined above contain more than two decimal places. If more 
decimal places are used, an error message will be printed. Execution may 
be continued but caution is advised, as invalid results may be generated. 

Columns 61: TYPE (l) - Label type (integer or real) to be 
placed on the x axis. A ’D' punched will plate integer numbers 
on the axis, and either an ’L' or ’C’ punched will place real 
numbers on the axis. 


Column 62: TYPE (2) - Same as TYPE (l) but for the y axis. 

Column 63: TYPE (3) - Type of plot to be generated (discrete, 
continuous, linear, or histogram). The options are specified 
as follows: 


Punch Option 

D This punch will cause the discrete option to be 

executed. Only the input data is plotted and the 
points are not connected in any manner. 

C This punch exercises the continuous plot option. 

After the individual points are plotted, the program 
calculates a best-fit function using the least- 
squares method and then plots the calculated 
function. (See Appendix A for details. ) 
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Punch 


Option 


L This punch causes the program to plot each indivi- 

dual point and at the same time connect each point 
to the previous one with a straight line. 

H This punch causes a histogram to be generated 

from the input. 

Column 64: NODAT - Punch any symbol in this column and no 
data points will be marked on the grid. 

• An example of this card is as follows: 

Columns 

1-10 11-20 21-30 31-40 41-50 51-60 61 62 63 64 

1.0 2.0 9.0 90.0 1.0 60.0 D C C * 

• Title Card 

This data card contains the labels for the grid. The format for 
this data card is (6A4, A2, 6A4, A2, 7A4). 

TITLX - This is a string of 26 characters to be placed on the 

x axis. 

TITLY - This is a string of 26 characters to be placed on the 

y axis. 

TITLE - This is a string of 28 characters to be placed above 
the grid. An example of the title card is as follows: 


Columns 

1-26 27-52 53-80 

PRESSURE C02 MM-HG EFFICIENCY-PCT. MOLECULAR SIEVE 


• Deck of x and y coordinates 

These data cards contain the x and y values of the data to be 
plotted and are punched according to the format 2F10.4. 
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When using the continuous or histogram option, the data points must 
be entered in the data set in order of increasing x value; that is, x must be 

less than x : x must be less than x . ... x must be less than x . If this 
2 2 3 n-1 n 

restriction is violated, invalid results will be obtained. 

• Data Card 


Columns 


1-10 11-20 

1.53 100.5 


• End Card 

An END card must be placed at the end of each data set. This 
card has the following form: 

Columns 

21-23 

END 

Multiple curves can be plotted on a single graph. Data for each new 
set can be stacked behind the previous case by adding a new CASE card and 
END card. Figure 1 illustrates a generalized GECAP deck setup. Section II 
shows a variety of different example plots which demonstrate the use of the 
different options available to the GECAP user. 


C. Operating Instructions 

The following is a list of operating instructions for using GECAP on the 
IBM 1130. 

1. Place disc in position in the disc storage unit. 

2. Turn file ON/OFF switch to the ON position. 

3. Wait for the FILE READY light to appear on the console typewriter. 

4. Press the NPRO button on card reader. 

5. Load the GECAP deck. 
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GENERALIZED GECAP DECK SET-UP 


000000000111111 1 11122222222223333333333444444444455555555556666666666777777 
12 34567890 12 3456 73901 2 34567890 12 345 67 890 12 343 6 7890123456 7 8901 23456789012345 
// JOB OR //***COLD START CARD***// 

// XEU GECAP 1 
•LOCALGECAP * LABEL*HIST »ERRO 
CASE CARD 1 

PROGRAM INSTRUCTION CARD 
TITLE CARD 

XVALUE ( 1 ) YVALUE ( 1 ) (THIS CARD BEGINS THE DATA SET FOR CASE 1) 

XVALUE ( 2 ) YVALUE ( 2 ) 


XVALUE ( N ) 

CASE CARD 
XVALUE! 1 ) 
XVALUE ( 2 ) 


XVALUE ( N ) YVALUE ( N ) 

END 


YVALUE ( N ) 

END 

2 

YVALUE ( 1 ) (THIS CARD BEGINS THE DATA SET FOR CASE 2) 

YVALUE (2 ) 


CASE CARD 
XVALUE ( 1 ) 
XVALUE ( 2 ) 


XVALUE t N I YVALUE ( N ) 

END 


—NOTE — MAXIMUM VALUE OF N IS 100. 

—NOTE — MAXIMUM VALUE OF M (NUMBER OF CASES) IS UNLIMITED. 

the number of data mark symbols, however* is 

LIMITED TO SIX. THEREFORE. IF DATA MARKS ARE USED. 
ONLY THE DATA POINTS FOR THE FIRST SIX CASES WILL 
BE LABELED WITH THESE MARKS. 

Figure 1. Generalized GECAP deck set-up. 


M 

YVALUE ( 1 ) (THIS CARD BEGINS THE DATA SET FOR CASE M) 

YVALUE ( 2 ) 


6. Check to see that all 15 sense-switches above the console type- 
writer are in the OFF position. 


7. Press the START button on card reader. 

8. If a Cold Start Card is being used, press the PROGRAM STOP 
button, press the RESET button, and then press PROGRAM LOAD button on the 
console typewriter. 


7 


9. If a JOB card is being used, press the green START button on 
the console typewriter. 

10. Wait for user message (l) to appear on the console typewriter. 
Press the START button on the console typewriter after this message has 
been written on the typewriter. 

11. Let the program run. ( Error messages may appear on the 
console typewriter indicating input errors) . 


SECTION III. EXAMPLE CASES 

Five example cases are shown in Figures 2 through 16 in an effort to 
illustrate the options and combinations of options available in GECAP. Each 
example includes the following information: 

1 . A listing of the input deck. 

2. A listing of the printed output. 

3. A copy of the output plot. 

Example Case 1 illustrates the linear option with data points marked. 
The x-axis is labeled with real numbers and the y-axis is labeled with integer 
numbers. Each of the plotted coordinates are connected with a straight line. 
Note that the output listing validates the value of the input variables. 


000000000 1111111111222222222233 i i i'i'ii 3 3 4 44 444 4 4 <♦ 4 b 335833b>3S66t>btoooOb 17/ II I I / I I a 
123456789012 34567890 12 34567890 1234S67S90 123456 7890123456 789012 34367890 12 34567890 
// JOB OR //***C0LD START CARO***// 

// XEU GECAP 1 
•LOCALGFCAP. LABEL •HIST tERRO 
CASE 1 OF X 

X. 5. 10. 100 . 0. 0 . COl 

TEST CASE X X-COORDINATE TEST CASE X Y-COORDINATE LINEAR OPTION 


0. 

0. 

1. 

XO. 

2. 

30. 

3. 

5. 

4. 

13.87 

5. 

81. 

6. 

XOO. 

7. 

54 . 

6 • 

49. 

9. 

49. 

XO. 

0.0 


END 


Figure 2. Listing of input deck for Example Case 1. 
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Example Case 2 illustrates the discrete option. Both axes are labeled 
with integer numbers. 

Example Case 3 illustrates the histogram option. Note that the maxi- 
mum x value for the case is specified as 20 on the Program Instruction Card 
while the maximum x value in the input data is only 19. When the H option 
is being used, the value for XMAX on the Program Instruction Card should 
always be at least one x increment larger than the maximum value of the 
input data. 

Example Case 4 illustrates the continuous option with no data points 
marked. Also illustrated here is the use of the exponential option to scale 
the axis values. Note that the x and y values of the input data are not scaled 
down but are within the specified limits of the program instruction card. 


Example Case 5 illustrates the continuous option with three cases 
plotted on the same graph. The plot includes symobls identifying each of the 
three corresponding curves. Notice that in the printed output, a user’s 
message appears which indicates that the polynomial calculated may not be 
accurate. After inspection of the curve on the output plot, it is seen that the 
curve is acceptable for most practical purposes. A printout of the polynomial 
coefficients is given for each curve along with the square of the errors from 
each data point. 


X-AXIS INCREMENT * 1.00 

X-AXIS LIMIT = 10.00 

X-AXIS ORIGIN = 0.00 

NCASES - 1 

TYPE! 1) » C 

TYPE t 2 ) * 0 

TYPE ( 3) » L 


Y-AXIS INCREMENT = 5.00 
Y-AXIS LIMIT = 100.00 
Y-AXIS ORIGIN * 0.00 


*** I NPUT DATA*** 

TEST CASE 1 X-C00RDINATE TEST CASE 1 Y-COORDINATE 


CASE 1 OF 1 


0.0000 

0.0000 

1.0000 

10.0000 

2.0000 

30.0000 

3.0000 

5.0000 

A. 0000 

13.8700 

5.0000 

81.0000 

6.0000 

100.0000 

7.0000 

5 A . 0000 

8.0000 

A9.0000 

9.0000 

A9.0000 

10.0000 

0.0000 


Figure 3. Listing of printed output for Example Case 1. 
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LINEAR OPTION 


8 



oivNraaoao-A t osvo isai 
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TE5T CASE 1 X- COORDINATE 

Figure 4. Plotted output for Example Case 1. 


00000000011111 111112222222222 J33333333:j44444444445355555355666666ij6o 67 77 77 77 7 / 7o 
1 2 34 567890 1 2 3456 7890 1 2 34 56 7 89P12 345 6 789 0 12 345 6 78 90 12 345 67 890 12 343 6 7 890 12 34 56 7d 90 
// JOB OR //***COLD START CARD***// 

// XEO GECAP 1 

•localgecap, label .hi st.erro 

CASE 1 OF 1 

2. 10. 20. 200. 0. 50. ODD 

TEST CASE 2 X-COORD I-NATE TEST CASE 2 Y-COORDINATE DESCRETE OPTION 

2 . 110 . 

4. 120. 

6. 115.56 

8. 170.2 

9,5 120. 

10. 125.5 

12. 130. 

14. 105. 

16. 190. 

18. 200. 

20. 160. 

END 

Figure 5. Listing of input deck for Example Case 2. 


X-AXIS INCREMENT ■ 

2*00 

Y-AXIS 

INCREMENT ■ 

10*00 

X-AXIS LIMIT ■ 

20*00 

Y-AXIS 

LIMIT ■ 

200*00 

X-AXIS ORIGIN • 

0*00 

Y-AXIS 

ORIGIN > 

50*00 

N CASES 

1 




TYPE ( 1) - 

D 




TYPE ( 2) ■ 

D 




TYPE! 3) ■ 

D 





♦♦♦INPUT DATA*** 

TEST CASE 2 X-COORDINATE TEST CASE 2 Y-COORDINATE 


CASE 1 OF 1 


2*0000 

4*0000 

6*0000 

8*0000 

9*5000 

10*0000 

12*0000 

14*0000 

16*0000 

18*0000 

20*0000 


110*0000 

120*0000 

115*5600 

170*2000 

120.0000 

125*5000 

130*0000 

105*0000 

190*0000 

200*0000 

160*0000 


Figure 6. Listing of printed putput for Example Case 2. 
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UOUOUUWIunaUl lli2222222222 3ii3.3333 33A‘.<.AA*.M**A‘o:35:>D5 5 33 votaoooooooo 7////////7 
123AS67B9012 3A56789012 3A5e789012 3A3678MO. l 23A56/B9ei.;JA5(>/o90i2J‘»50i'>,'»ui.3‘OU/oy 
// jos or / / ***cold sr*t>r card***// 

// XEU GECAP 1 
•LOCAL6FCAP.LABFL.HIST .EH WO 
CASK 1 OF 1 

1. 20 . 20# i?J. im £.J« L/Clt 

TEST CASE 3 — CELL MEAN Tt.iT C Ajt 3 — FHEOUEnLY Hl.iTJC.WAV uPIlu.Ti 

1 . 25. 

2 . 39 . 

3. 50. 

a . 55 . 

5. 58. 

6 • 6 0 . 

7. 75. 

8. 95. 

9. 103.98 

10. 115. 2A 

11. 116.5 

12 . 101 . 

13. 92. 

1 A » 8 A , 

15. 60. 

16. A2.25 

17. 36. 

18. 30. 

19. 22.0A 

END 

Figure 8. Listing of input deck for Example Case 3. 


X-AXIS INCREMENT ■ 

1.00 

Y-AXIS 

INCREMENT * 

20.00 

X-AXIS LIMIT = 

20.00 

Y-AXIS 

LIMIT = 

120.00 

X-AXIS ORIGIN * 

1.00 

Y-AXIS 

ORIGIN = 

20.00 

NCASES 

1 




TYPE! 1) * 

D 




TYPE! 2) * 

C 




TYPE! 3) = 

H 




•••INPUT DATA*** 





TEST CASE 3— CELL MEAN TEST CASE 3— FREOUENCY 


CASE 1 OF 1 


1.0000 

2.0000 

3.0000 
A. 0000 

5.0000 

6.0000 

7.0000 

8.0000 
9.0000 

10.0000 

11.0000 

12.0000 

13.0000 
1A.0000 

15.0000 

16.0000 

17.0000 

18.0000 
19.0000 


25.0000 

39.0000 

50.0000 

55.0000 

58.0000 

60.0000 

75.0000 

95.0000 
103.9800 
115.2A00 
116.5000 
101.0000 

92.0000 
8 A • 0000 

60.0000 
A2.2500 

36.0000 

30.0000 
22.0A00 



Figure 9. Listing of printed output for Example Case 3. 
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HISTOGRAM OPTION 



Sj 8 8 8 ? 8 

AGN3nD3dJ--E 3SV0 1S31 
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TEST CASE 3- -CELL MEAN 

Figure 10. Plotted output for Example Case 3. 




OOOOOOOOOlllli 11 1 112222222222333J iiiii 3 4‘.4'.<.<.4^‘,i53353533a6bl)t)o6ofa6t) 7/7/77 / '///b 
123 4567 890 12 34 56 7890 12345678901 2 34b 6 789 01 2 3456 7890 l 2 34 5 6 789012 34 56 7 6901234567890 
// JOB OR //***COLO START CARD»**// 

// XEU GECAP 1 

*LOC ALGEC AP • L ARE L » H I S T tERRO 

CASE 1 OF 1 -1 -1 

5. 5. 40. 10. 0. -10. CCC* 

TEST CASE 4 X-COORD I NA TE TEST CASE 4 Y-COORD I (MAT E CONT. OPTION. NU DATA MARKS 

0 . 10 . 

1.05 9.95 

I. 92 9.82 

2.79 9.61 

3.67 9.34 

4.54 8.99 

5.4] 8.57 

6.28 8.09 

7.16 7.55 

8.03 6.95 

8.90 6.29 

9.77 5.59 

10.65 4.85 

II. 52 4.07 

12.39 3.26 

13.26 2.42 

14.14 1.56 

15.01 0.70 

15.71 0.0 

16.76 - 1.05 

17.63 - 1.91 

18.50 - 2.76 

19.38 - 3.58 

20.25 - 4.38 

21.12 - 5.15 

21.99 - 5.88 

22.86 - 6.56 

23.73 - 7.19 

24.60 - 7.77 

25.47 - 8.29 

26.34 - 8.75 

27.21 - 9.14 

28.08 - 9.46 

28.85 - 9.70 

29.72 - 9.88 

30.59 - 9.98 

31.46 - 10.00 

END 


Figure 11. Listing of input deck for Example Case 4. 
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X-AXIS INCREMENT * 


o 

o 

• 

m 

Y-AX 

IS 

X-AXIS LIMIT « 


40.00 

Y-AX 

IS 

X-AXIS ORIGIN ■ 


0.00 

Y-AX 

IS 

NCASES • 


1 



TYPE ( 1) • 

C 




TYPE ( 2) - 

C 




TYPE I 3) ■ 

C 





*** I NPUT 

data*** 

TEST CASE 4 X-COORDINATE 

TEST CASE 


CASE 

1 OF 1 

0.0000 


10.0000 

1.0500 


9.9500 

1.9200 


9.8200 

2.7900 


9.6100 

3.6700 


9.3400 

4.5400 


8.9900 

5.4100 


8.5700 

6.2800 


8.0900 

7.1600 


7.5500 

8.0300 


6.9500 

8.9000 


6.2900 

9.7700 


5.5900 

10.6500 


4.8500 

11.5200 


‘ 4.0700 

12.3900 


3.2600 

13.2600 


2.4200 

14.1400 


1.5600 

15.0100 


0.7000 

15.7100 


0.0000 

16.7600 


-1.0500 

17.6300 


-1.9100 

18.5000 


-2.7600 

19.3800 


-3.5800 

20.2500 


-4.3800 

21.1200 


-5.1500 

21.9900 


-5.8800 

22.8600 


-6.5600 

23.7300 


-7.1900 

24.6000 


-7.7700 

25.4700 


-8.2900 

26.3400 


-8.7500 

27.2100 


-9.1400 

28.0800 


-9.4600 

28.8500 


-9.7000 

29.7200 


-9.8800 

30.5900 


-9.9800 

31.4600 


-10.0000 


INCREMENT - 
LIMIT * 
ORIGIN - 


4 Y-COORDINATE 


Figure 12. Listing of printed output for Example Case 4 
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•CO-EFFICIENTS FOR POLYNOMIAL OF DEGREE 
10.000709097832 
0.003862024605 
—0.051973071545 
0.000391693861 
0.000002719892 
0. 000002066409 
-0.000000070663 
0.000000000697 


•CALCULATED VALUES AND SUM OF ERRORS SQUARED FOR 7 DEGREE POLYNOMIAL* 


FOR 

POINT 

1 

YCALC 

= 

10.000709 

D I FF**2 

Z 

0.000000502803 

FOR 

POINT 

2 

YCALC 

S 

9.947923 

D I FF**2 

= 

0.000004313210 

FOR 

POINT 

3 

YCALC 

= 

9.819390 

DIFF**2 

s 

0.000000371601 

FOR 

POINT 

4 

YCALC 

s 

9.615908 

D1FF**2 

= 

0.000034915439 

FOR 

POINT 

5 

YCALC 

= 

9.335927 

DIFF**2 

= 

0.000016588005 

FOR 

POINT 

6 

YCALC 

x 

8.988198 

D I FF**2 

= 

0.000003247140 

FOR 

POINT 

7 

YCALC 

s 

8.572699 

DIFF**2 

X 

0.000007289361 

FOR 

POINT 

8 

YCALC 

x 

8.092588 

D I FF **2 

3 

0.000006702353 

FOR 

POINT 

9 

YCALC 

s 

7.544891 

DIFF**2 

S 

0.000026091918 

FOR 

POINT 

10 

YCALC 

X 

6.946120 

D I FF**2 

= 

0.000015053024 

FOR 

POINT 

11 

YCALC 

s 

6.294849 

D I FF**2 

3 

0.000023513143 

FOR 

POINT 

12 

YCALC 

s 

5.595945 

D 1 FF**2 

S 

0.000035346223 

FOR 

POINT 

13 

YCALC 

s 

4.845912 

D I FF**2 

= 

0.000016704003 

FOR 

POINT 

14 

YCALC 

s 

4.067449 

D I FF**2 

= 

0.000006506254 

FOR 

POINT 

15 

YCALC 

= 

3.258165 

D I FF **2 

3 

0.000003365522 

FOR 

POINT 

16 

YCALC 

X 

2.424215 

D I FF**2 

3 

0.000017769932 

FOR 

POINT 

17 

YCALC 

3 

1.562073 

D I FF**2 

S 

0.000004300971 

FOR 

POINT 

18 

YCALC 


0.697891 

D I FF ** 2 

- 

0.000004445361 

FOR 

POINT 

19 

YCALC 

- 

-0.001492 

D IFF**2 

* 

0.000002226602 

FOR 

POINT 

20 

YCALC 

B 

-1.049323 

D I FF**2 

3 

0.000000457832 

FOR 

POINT 

21 

YCALC 

B 

-1.909124 

D I FF**2 

X 

0.000000767067 

FOR 

POINT 

22 

YCALC 

X 

-2.754392 

DI FF **2 

3 

0.000031449202 

FOR 

POINT 

23 

YCALC 

= 

-3.588137 

DIFF**2 

X 

0.000066218449 

FOR 

POINT 

24 

YCALC 

3 

-4.385229 

D I FF**2 

E 

0.000027344544 

FOR 

POINT 

25 

YCALC 

• 

-5.149389 

D I FF**2 

3 

0.000000372938 

FOR 

POINT 

26 

YCALC 

E 

-5.875004 

DIFF**2 

= 

0.000024950876 

FOR 

POINT 

27 

YCALC 

r 

-6.556761 

D I FF**2 

E 

0.000010490893 

FOR 

POINT 

28 

YCALC 

B 

-7.189656 

D I FF**2 

X 

0.000000118064 

FOR 

POINT 

29 

YCALC 

E 

-7.769007 

D I FF**2 

= 

0.000000985816 

FOR 

POINT 

30 

YCALC 

X 

-8.290437 

DIFF**2 

Z 

0.000000191150 

FOR 

POINT 

31 

YCALC 

3 

-8.749856 

DI FF**2 

X 

0.000000020600 

FOR 

POINT 

32 

YCALC 

X 

-9.143423 

D I FF**2 

3 

0.000011718138 

FOR 

POINT 

33 

YCALC 

X 

-9.467486 

D I FF**2 

X 

0.000056053369 

FOR 

POINT 

34 

YCALC 

E 

-9.693488 

D I FF**2 

3 

0.000042395776 

FOR 

POINT 

35 

YCALC 

X 

-9.876957 

DI FF **2 

3 

0.000009259824 

FOR 

POINT 

36 

YCALC 

3 

-9.980715 

D I FF**2 

X 

0.000000511452 

FOR 

POINT 

37 

YCALC 

X 

-10.000982 

D IFF**2 

= 

0.000000964853 


SUM OF SQUARES OF ERRORS FOR CURVE-FIT = 0.000513523724 

Figure 12. (Concluded) 
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OOOOOOOUO 1 11 li 1 1 i 1 12222J2^2 2 V 3 33 3 3 3 5 3 3 3<.<.4‘.h<«‘.<.<»4 ;'j5iini‘j-j36ot>ocfaotioo7 17 / / / II I la 
1234S67B901234567S9J1/3456 U23456 /i‘». 123456 7 '»s . 123450 74901 23456 7890 j.2 3456 7H9U 

// JOB OR //***COLD START CARD***// 

// XEO GECAP 1 
*LOCALGECAP. LABEL .HIST .ERRO 
CASE 1 OF 3 EXAMPLE CURVE 1 

.1 .10001 .8 1.1 0.0 o.o CCC 

TEST CASE 5 X-COORDINATE TEST CASE 5 Y -COORD I NA T b MULTIPLE CONTINUOUS UP I ION 
O.OAOOO 0.00164 
0.06000 0.00883 

0.08000 0.027S8 

0.10000 0.06336 

0.12000 0.11892 

0.14000 0.19338 

0.16000 0.28254 

0.18000 0.38022 

0.20000 0.47972 

0.22000 0.57508 

0.24000 0.66188 

0.26000 0.73748 

0.28000 0.80085 

0.30000 0.85216 

0.32000 0.89247 

0.34000 0.92331 

0.36000 0.94630 

0.38000 0.96305 

0.40000 0.97501 

0.42000 0.98337 

0.44000 0.98913 

0.46000 0.99301 

0.43000 0.99560 

END 

CASE 2 OF 3 EXAMPLE CURVE 2 
0.12000 0.01022 
0.14000 0.02879 

0.16000 0.05832 

0.18000 0.10040 

0.20000 0.15527 

0.22000 0.22173 

0.24000 0.29730 

0.26000 0.37858 

0.28000 0.46193 

0.30000 0.54363 

0.32000 0.62076 

0.34000 0.69098 

0.36000 0.75305 

0.38000 0.80622 

0.40000 0.85065 

0.42000 0.88680 

0.44000 0.91556 

0.46000 0.93806 

0.48000 0.95526 

0.50000 0.96817 

0.52000 0.97769 

0.54000 0.98463 

0.56000 0.98959 

0.58000 0.99311 

0.60000 0.99555 

END 


Figure 14. Listing of input deck for Example Case 5. 
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CASE 3 OF 3 
0.22000 0 
0.24000 0 

0.26000 0 
0.28000 0 
0.30000 0 

0.32000 0 

0.34000 0 

0.36000 0 

0.38000 0 

0.40000 0 

0.42000 0 

0.44000 0 

0.46000 0 

0.48000 0 

0.50000 0 

0.52000 0 

0.54000 0 

0.56000 0 

0.58000 0 

0.60000 0 
0.62000 0 
0.64000 0 

0.66000 0 
0.68000 0 
0.70000 0 

0.72000 0 

0.74000 0 

0.76000 0 

0.78000 0 

0.80000 0 


EXAMPLE CURVE 3 
01197 
03007 
05566 
08968 
13234 
18334 
24188 
30617 
37426 
4441 5 
51360 
58064 
64370 
70159 
75377 
79961 
83909 
87246 
90018 
92280 
94100 
95543 
96671 
97540 
98202 
98698 
99067 
99338 
99535 
99677 

END 


Figure 14. (Concluded). 
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X-AXIS INCREMENT ■ 


0*10 

Y-AXIS 

INCREMENT - 

X-AXIS LIMIT « 


0.80 

Y-AXIS 

LIMIT » 

X-AXIS ORIGIN - 


0.00 

Y-AXIS 

ORIGIN « 

NCASES 

9 


3 



TYPE ( 1) 

9 

C 




TYPE ( 2) 

U 

c 




TYPE < 3) 

m 

c 





♦♦♦INPUT 

data*** 


TEST CASE 

5 X-COORDINATE 

TEST CASE 

5 Y— COORO I NATE 


CASE 1 OF 3 


0.0400 

0.0016 

0.0600 

0.0088 

0.0800 

0.0275 

0.1000 

0.0633 

0.1200 

0.1189 

0. 1400 

0.1933 

0.1600 

0.2825 

0.1800 

0.3802 

0.2000 

0.4797 

0.2200 

0.5750 

0.2400 

0.6618 

0.2600 

0.7374 

0.2800 

0.8008 

0.3000 

0.8521 

0.3200 

0.8924 

0.3400 

0.9233 

0.3600 

0.9463 

0.3800 

0.9630 

0.4000 

0.9750 

0.4200 

0.9833 

0.4400 

0.9891 

0.4600 

0.9930 

0.4800 

0.9956 

Figure 15. 

Listing of printed output for 1 


O.IO 

1*10 

0*00 
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♦CO-EFFICIENTS FOR POLYNOMIAL OF DEGREE 9* 
-0.039928562066 
2.680569725111 
-65.717747807502 
718.222608039447 
—2499 .070896148681 
1838.338735103607 
6513.944360733032 
-10910.395404815673 
-2377.100768089294 
9138.639202117919 


♦CALCULATED VALUES AND SUM OF ERRORS SQUARED FOR 9 DEGREE POLYNOMIAL* 


FOR 

POINT 

1 

YCALC 

x 

0.001927 

D I FF **2 

= 

G.U00000082703 

FOR 

POINT 

2 

YCALC 

a 

0.008772 

D I FF**2 

X 

0.000000003316 

FOR 

POINT 

3 

YCALC 

a 

0.026791 

D I FF**2 

X 

0.000000621963 

FOR 

POINT 

4 

YCALC 

s 

0.063058 

D I FF **2 

a 

0.000000091133 

FOR 

POINT 

5 

YCALC 

= 

0.119515 

D I FF **2 

x 

0.000000354491 

FOR 

POINT 

6 

YCALC 

= 

0.194297 

D I FF **2 

X 

0.000000841481 

FOR 

POINT 

7 

YCALC 

X 

0.283006 

D I FF**2 

X 

0.000000217405 

FOR 

POINT 

8 

YCALC 

a 

0.379886 

D I FF**2 

X 

0.000000111237 

FOR 

POINT 

9 

YCALC 

s 

0.478844 

D I FF**2 

X 

0.000000766828 

FOR 

POINT 

10 

YCALC 

= 

0.574275 

DIFF**2 

= 

0.000000647926 

FOR 

POINT 

11 

YCALC 

= 

0.661666 

D I FF**2 

X 

0.000000045699 

FOR 

POINT 

12 

YCALC 

X 

U. 737961 

D I FF**2 

X 

0.000000231697 

FOR 

POINT 

13 

YCALC 

s 

0.801681 

D I FF**2 

a 

0.000000691816 

FOR 

POINT 

14 

YCALC 

s 

0.852820 

D 1 FF**2 

X 

0.000000436204 

FOR 

POINT 

15 

YCALC 

s 

0.892541 

D I FF**2 

X 

0.000000005145 

FOR 

POINT 

16 

YCALC 

s 

0.922735 

DIFF**2 

a 

0.000000330597 

FOR 

POINT 

17 

YCALC 

X 

0.945499 

D IFF**2 

X 

0.000000640990 

FOR 

POINT 

18 

YCALC 

s 

0.962653 

D 1 FF**2 

X 

0.000000157419 

FOR 

POINT 

19 

YCALC 

X 

0.975387 

D I FF**2 

3 

0.000000142702 

FOR 

POINT 

20 

YCALC 

X 

0.984216 

D I FF **2 

X 

0.000000716604 

FOR 

POINT 

21 

YCALC 

3 

0.989394 

DIFF**2 

3 

0.000000069805 

FOR 

POINT 

22 

YCALC 

3 

0.992002 

D I FF**2 

3 

0.000001015472 

FOR 

POINT 

23 

YCALC 

= 

0.995954 

D 1 FF**2 

= 

0.000000125563 

SUM 

OF SQUARES 

i OF ERRORS 

FOR CURVE-FIT « 

0.000008348206 



Figure 15. (Continued). 
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*** I NPU T DATA*** 


TEST CASE 


***USERS 


5 X-COORDINATE TEST CASE 5 Y— COORD 1 NATE 


CASE 2 OF 3 


O. 1200 

0.0102 

0. 1 LOO 

0.0287 

0.1600 

0.0583 

0.1800 

0.1004 

0.2000 

0.1652 

0.2200 

0.2217 

0. 2400 

0.2973 

0.2600 

0.3785 

0.2800 

0.4619 

0. 3000 

0.5436 

0.3200 

0.6207 

0.3400 

0.6909 

0.3600 

0.7530 

0.3800 

0.8062 

0.4000 

0.8506 

0.4200 

0.8868 

0.4400 

0.9155 

0.4600 

0.9380 

0.4800 

0.9552 

0.6000 

0.9681 

0.6200 

0.9776 

0.6400 

0.9846 

0.6600 

0.9895 

0.6800 

0.9931 

0.6000 

0.9955 


MESSAGE*** THE POLYNOM I NAL CALCULATED MAY OR MAY NOT BE THE BEST FITTED 
CURVE 


Figure 15. (Continued). 
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•CO-EFFICIENTS FOR POLYNOMIAL OF DEGREE 9* 
0.733166191959 
-17.548953562974 
157.079668283462 
-618.872270822525 
296.712045907974 
9378.129894256591 
-44079.530563354492 
90706.581665039062 
-90827.897857666015 
36035.444580078125 


•CALCULATED VALUES AND SUM OF ERRORS SQUARED FOR 9 DEGREE POLYNOMIAL* 


FOR 

POINT 

1 

YCALC 

• 

0.011873 

D I FF**2 

= 

0.000002733528 

FOR 

POINT 

2 

YCALC 

M 

0.026310 

D I FF**2 

= 

0.000006146741 

FOR 

POINT 

3 

YCALC 

U 

0.056916 

D 1 FF **2 

a 

0.000001969616 

FOR 

POINT 

4 

YCALC 

X 

0.101147 

D I FF**2 

a 

0.000000558885 

FOR 

POINT 

5 

YCALC 

■ 

0. 157209 

D I FF**2 

= 

0.000003763071 

FOR 

POINT 

6 

YCALC 

a 

0.223401 

D I FF**2 

a 

0.000002793994 

FOR 

POINT 

7 

YCALC 

a 

0.297697 

D I FF **2 

a 

0.000000157887 

FOR 

POINT 

8 

YCALC 

a 

0.377576 

D I FF**2 

a 

0.000001006416 

FOR 

POINT 

9 

YCALC 

a 

0.460065 

D I FF**2 

a 

0.000003477294 

FOR 

POINT 

10 

YCALC 

m 

0.541940 

DIFF*«2 

a 

0.000002853150 

FOR 

POINT 

11 

YCALC 

Z 

0.620046 

DIFF**2 

a 

0.000000509448 

FOR 

POINT 

12 

YCALC 

r 

0.691631 

D I FF**2 

X 

0.000000425034 

FOR 

POINT 

13 

YCALC 

a 

0.754671 

D I FF **2 

a 

0.000002630334 

FOR 

POINT 

14 

YCALC 

a 

0.808083 

D IFF **2 

a 

0.000003472011 

FOR 

POINT 

15 

YCALC 

a 

0.861795 

D IFF**2 

a 

0.000001312294 

FOR 

POINT 

16 

YCALC 

a 

0.886659 

D I FF**2 

a 

0.000000019819 

FOR 

POINT 

17 

YCALC 

a 

0.914176 

D I FF**2 

a 

0.000001914500 

FOR 

POINT 

18 

YCALC 

a 

0.936084 

D I FF**2 

a 

0.000003901386 

FOR 

POINT 

19 

YCALC 

a 

0.953881 

D I FF**2 

a 

0.000001901406 

FOR 

POINT 

20 

YCALC 

a 

0.968385 

D I FF**2 

a 

0.000000046550 

FOR 

POINT 

21 

YCALC 

a 

0.979562 

D I FF**2 

a 

0.000003508122 

FOR 

POINT 

22 

YCALC 

s 

0.986698 

D I FF**2 

a 

0.000004279220 

FOR 

POINT 

23 

YCALC 

a 

0.989476 

D I FF**2 

a 

0.000000012949 

FOR 

POINT 

24 

YCALC 

z 

0.990075 

D IFF **2 

a 

0.000009207652 

FOR 

POINT 

25 

YCALC 

■ 

0.996855 

D I FF**2 

a 

0.000001704946 


SUM OF SQUARES OF ERRORS FOR CURVE-FIT = 0.000060306266 

Figure 15. ( Continued) . 
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*** I NPUT DATA*** 


CASE 5 X-COORDINATE 

TEST CASE 


CASE 3 OF 

3 

0.2200 


0.0119 

0.2400 


0.0300 

0.2600 


0.0556 

0.2800 


0.0896 

0. 3000 


0.1323 

0.3200 


0.1833 

0.3400 


0.2418 

0.3600 


0.3061 

0.3800 


0.3742 

0.4000 


0.4441 

0.4200 


0.5136 

0.4400 


0.5806 

0.4600 


0.6437 

0.4800 


0.7015 

0. 5000 


0.7537 

0.5200 


0.7996 

0.5400 


0.8390 

0. 5600 


0.8724 

0.5800 


0.9001 

0.6000 


0.9228 

0.6200 


0.9410 

0.6400 


0.9554 

0.6600 


0.9667 

0.6800 


0.9754 

0.7000 


0.9820 

0.7200 


0.9869 

0.7400 


0.9906 

0.7600 


0.9933 

0.7800 


0.9953 

0.8000 


0.9967 


Figure 15. (Continued). 


Y-COORDINATE 
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•CO-EFFICIENTS FOR POLYNOMIAL OF DEGREE 9* 

-0.848485231399 
20.262203224003 
-197.173971712589 
10 59. 8804 A 1665649 
-3750.190005111694 
9700.356727600097 
-17520.517059326171 
20004.520645141601 
-12707.314575195312 

3398.905395507812 » 


•CALCULATED 

VALUES 

AND SUM OF 

ERRORS SQUARED 

FOR 9 DEGREE 

POLYNOMIAL* 

FOR 

POINT 

1 

YCALC 

Z 

0.012745 

D I FF **2 

= 

0.000000601850 

FOR 

POINT 

2 

YCALC 

z 

0.029356 

D I FF**2 

a 

0.000000509098 

FOR 

POINT 

3 

YCALC 

= 

0.054756 

DIFF**2 

3 

0.000000815690 

FOR 

POINT 

4 

YCALC 

z 

0.089296 

DIFF**2 

z 

0.000000146885 

FOR 

POINT 

5 

YCALC 

* 

0. 132728 

DIFF**2 

z 

0.000000151046 

FOR 

POINT 

6 

YCALC 

z 

0.184288 

D I FF**2 

z 

0.000000899720 

FOR 

POINT 

7 

YCALC 

z 

0.242783 

DIFF**2 

z 

0.000000816801 

FOR 

POINT 

8 

YCALC 

z 

0.306682 

D I FF**2 

z 

0.000000262453 

FOR 

POINT 

9 

YCALC 

z 

0.374215 

D I FF**2 

s 

0.000000002024 

FOR 

POINT 

10 

YCALC 

z 

0.443483 

DIFF**2 

z 

0.000000444748 

FOR 

POINT 

11 

YCALC 

z 

0.512575 

D I FF**2 

= 

0.00000 10504O7 

FOR 

POINT 

12 

YCALC 

z 

0.579674 

DIFF**2 

8 

0.000000933148 

FOR 

POINT 

13 

YCALC 

c 

0.643161 

D I FF**2 

a 

0.000000289499 

FOR 

POINT 

14 

YCALC 

B 

0.701712 

D IFF **2 

a 

0.000000015045 

FOR 

POINT 

15 

YCALC 

z 

0.754354 

D I FF**2 

a 

0.000000341927 

FOR 

POINT 

16 

YCALC 

z 

0.800516 

DIFF**2 

a 

0.000000821663 

FOR 

POINT 

17 

YCALC 

z 

0.840033 

DIFF**2 

a 

0.000000889492 

FOR 

POINT 

18 

YCALC 

z 

0.873110 

D I FF**2 

z 

0.000000423290 

FOR 

POINT 

19 

YCALC 

z 

0.900291 

DIFF**2 

a 

0.000000012349 

FOR 

POINT 

20 

YCALC 

B 

0.922329 

D I FF*»2 

m 

0.000000221075 

FOR 

POINT 

21 

YCALC 

z 

0.940093 

D I FF**2 

a 

0.000000822142 

FOR 

POINT 

22 

YCALC 

a 

0.954434 

DIFF**2 

B 

0.000000990398 

FOR 

POINT 

23 

YCALC 

a 

0.966066 

DIFF**2 

z 

0.000000414041 

FOR 

POINT 

24 

YCALC 

z 

0.975456 

DIFF**2 

z 

0.000000003223 

FOR 

POINT 

25 

YCALC 

a 

0.982819 

DIFF**2 

3 

0.000000639291 

FOR 

POINT 

26 

YCALC 

a 

0.988129 

D IFF**2 

Z 

0.000001321069 

FOR 

POINT 

27 

YCALC 

z 

0.991351 

DIFF**2 

S 

0.000000464910 

FOR 

POINT 

28 

YCALC 

B 

0.992795 

DJFF*«2 

Z 

0.000000341957 

FOR 

POINT 

29 

YCALC 

a 

0.993799 

DIFF**2 

8 

0.000002402732 

FOR 

POINT 

30 

YCALC 

3 

0.997637 

DIFF**2 

Z 

0.000000752988 

SUM 

OF SQUARES OF ERRORS FOR 

CURVE-FIT - 

0.000017800973 



Figure 15. (Concluded). 
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MULTIPLE CONTINUOUS OPTION 



31VNICMX): 


TEST CA5E 5 X- COORDINATE 

Figure 16. Plotted output for Example Case 5 



SECTION IV. PROGRAM MESSAGES 


A. User Messages 

The following two user-messages have been included in GECAP. 

1 »«** POSITION PLOTTER PEN APPROXIMATELY 3 INCHES FROM THF **** 

• • •** RIGHT EDGE OF THE PLOTTER. . . . PRESS START . . . **** 

This message appears on the console typewriter at the beginning of 
each execution of the program. Its purpose is to instruct the user to make 
sure the plotter pen is in a proper position relative to the plotting paper. To 
continue execution the green PROGRAM START button on the console type- 
writer must be pressed. (The plot origin is established at this position. ) 

2. ***USERS MESSAGE*** THE POLYNOMIAL CALCULATED MAY 
OR MAY NOT BE THE BEST FITTED CURVE 

This message occurs when a statistical test in BESFT has not been 
satisfied. Execution will continue and the user may determine if the calculated 
polynomial is acceptable for his purposes. (See the section on program limi- 
tations for more information. ) 


B. Error Messages 

GECAP was designed such that input to the program was made as simple 
as possible. However, various restrictions within the program logic may 
result in erroneous output if invalid input is used. In order to aid the GECAP 
user in detecting these errors, various messages are printed on the console 
typewriter, which indicate some of the common errors that occur from 
incorrect program instructions or input data. The following list describes 
these messages and gives the course of action taken by GECAP during their 
occurrence. 

1. ••••ERROR**** CASE CARD INPUT VARIABLE I* I NCORRFCT 

•NUMfER OF CASES HAS BEEN SET EQUAL TO 1 ^EXECUTION RESUMED* 
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This message indicates that an error has been made in the format of 
the CASE card. Execution will continue, setting the number of plots equal 
to 1. 


**“ERROR**“END POINT OF X AXIS IS LESS THAN STARTING POINT**** 
•“‘EXECUTION DISCONTINUED**** 

This message indicates that an error has been made when specifying 
the limits for the x axis. Execution will be discontinued at this point. 

3 * **** ERROR**** END POINT OF Y AXIS IS LESS THAN STARTING POINT**** 

““EXECUTION DISCONTINUED**** 

This message indicates the same thing as message 4 except for the 
y axis. Execution will be discontinued. 

4. ***ERROR“* NUMBER OF TIC MARKS EXCEEDS LEGIBILITY LIMIT*** 
STANDARD FIX-UP TAKEN, EXECUTION CONTINUING 

This message indicates that the values specified for labeling the axes 
calls for more than 20 ’tic’ marks or labeling intervals. The program will 
internally readjust the increments so that 10 'tic' marks are placed on the 
axes and scale values are recalculated. Execution will continue. The user 
may increase the increment value XINC or YINC to eliminate this error. 

5. •••• ERROR. . MORE THAN 100 DATA POINTS ON INPUT **** 

**** EXECUTION CONTINUING WITHOUT REMAINING POINTS** 

A limit of 100 data points has been established for this program. 

If this restriction is violated, the program will discard the extra data points 
and continue execution with the first 100 values. This message will be printed 
out for each input data card over 100. 

****ERR0R****THE VALUES SPECIFIED FOR LABELING THE AXES 

REQUIRE MORE PRECISION THAN THAT ALLOWABLE IN GECAP 

...CHECK GECAP USER MANUAL FOR DETAILS... 

THIS ERROR MAY BE INSIGNIFICANT IN SPECIAL CASES.. .IF EXECUTION 
IS STILL DESIRED PRESS START. ..IF NOT, PRESS STOP 
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This message indicates that more than two decimal places were 
requested on the program instruction card. In some cases this error will 
cause invalid results. If execution is desired, the green START button on 
the console typewriter must be pressed. Six different input parameters 
(XINC, YINC, XMAX, YMAX, XSTRT, YSTRT) are checked by the program 

and any one of these values could cause this message to appear. The START * 

button must be pressed each time the message appears, and, only after all 
six values have been checked, will execution of GECAP be resumed. 

7. ****ERROR****A data point was found to exceed the limits of the axes...**** 

♦•♦♦EXECUTION DISCONTINUED**** 


This message indicates that one of the input coordinates does not 
fall within the user defined limits of the coordinate axes. Execution will be 
discontinued. 


**** ERROR**** INCORRECT PLOT OPTION WAS USED FOR 

VARIABLE TYPE(l) .*♦** CHECK GECAP USER MANUAL FOR DFTAILS 

****TYPE IN THE DESIRED OPTION FROM THE CONSOLE TYPEWRITER 
(C,D,H, OR L). ****PRESS (EOF) BUTTON TO CONTINUE EXECUTION***** 


9. ****ERR0R**** INCORRECT PLOT OPTION WAS USED FOR 

VARIABLE TYPE(2).**** CHECK GFCAP USER MANUAL FOR DFTAILS. 
****TYPE IN THE DESIRED OPTION FROM THE CONSOLE TYPEWRITER 
(C,D,H, OR L). **** PRESS (EOF) BUTTON TO CONTINUE EXECUTION 


10. ****ERR0R**** INCORRECT PLOT OPTION WAS USED FOR 

VARIABLE TYPE(3) .***♦ CHECK GECAP USER MANUAL FOR DFTAILS. 
****TYPE IN THE DESIRED OPTION FROM THE CONSOLE TYPEWRITER 
(C,D,H, OR L). ****PRESS (EOF) BUTTON TO CONTINUE EXECUTION 


Messages 8, 9, and 10 appear when incorrect symbol options were 
placed in columns 61, 62, or 63 of the Program Instruction Card. When this 
message appears, the correct symbol should be typed in from the console 
typewriter. To resume execution, the EOF button on the typewriter should be 
pressed. 
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***EKROR*** ERROR IN ROUTINE CURF. EXECUTION CONTINUING 


This error message occurs in routine CURF and indicates that an 
error resulted while calculating the coefficients for the least squares poly- 
nominal. Execution will continue. 


SECTION V. PROGRAM DESCRIPTION 
A. General 


GECAP was designed exclusively for the IBM 1130 computing system. 
This system was chosen because it is common among industrial and govern- 
mental facilities, and therefore, provides easy access to the user. Core 
requirements for this program are 8K words of storage. 5 

This program has been constructed to provide the user with a rapid 
and accurate method of plotting a set of data on a grid without having to write 
any programming instructions. The program reads in information concerning 
starting values for the coordinate axes, upper limits for the axes, and step 
increments between ’TIC’ marks for the axes. These values are used to 
calculate scale values for the grid. The scale values calculated by the pro- 
gram are based on the assumption that the plot is small enough to fit an 
8. 5 x 11 in. sheet of paper. User titles are read in for the x and y axes and 
also for the grid, and placed at the proper location on the graph. 

Several different types of plots are available. The user selects the 
type of plot, the program generates the requested graph for the given data 
set. The supplied information and the input data is written on the printer for 
user reference. 

Numerous options concerning the type of plot and the method of labeling 
the axes are available to the user. 

Labeling the axes: 

1. Either or both of the coordinate axes may be labeled with integer 
or real numbers. 


5. The program is actually larger than 8K but by the use of the LOCAL and 
LINK system overlay routines, the program is within the 8K limitation of the 
1130 computer system. 
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2. The user may wish to indicate that his labeled values on the 
axes are scaled to a power of 10. If so, an option is available which will 
label the last 'TIC’ mark on each axis in such a manner as to indicate 
this scaling factor (ex: 20 x 10 _1 ). 

Type of plot: 

1. The data may be plotted as individual points marked by symbols 
on the grid. 

2. The data may be plotted as individual points with straight lines 
connecting each point to the previous point. 

3. A histogram option is available which generates a histogram 
from the user-supplied frequency information. 

4. An option may be used which plots the points individually and 
then performs a least-squares curve fit on the supplied data. The program 
uses a statistical test 6 to determine which degree polynomial best fits the 
data and then plots the curve. The coefficients of the calculated polynomial 
along with the errors associated with the differences from the calculated 
curve and the actual data are displayed on the printout. 

5. The user may not desire data point symbols on his graph. If so, 
an option is available, applicable to any of the above options, which causes no 
data point marks to be placed on the grid. 

6. Multiple cases may be plotted on the same grid for options 1, 2, 
4, and 5. 

7. If data point marks are used, one permits the data symbols to be 

labeled. 

After each set of data is plotted on the grid, the plotter pen is moved 
to a point below the x axis. The symbol associated with that set of data is 
drawn and followed by a user-specified name for that particular data set. 

Only six curves may be labeled in this fashion as the IBM 1130 plotter routines 
are limited to six different symbols for data point marks. 


6. See Appendix A . 
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B. Program Structure 


GECAP is built primarily in two sections. The mainline, GECAP, and 
supporting routines LABL, ERRO, and HIST perform the operations concerned 
with drawing and labeling the grid, plotting the individual data points, and 
generating a plot for all options except the least-squares curve-fit option. 

The curve-fit option is contained in the second main routine, BESFT, which 
calculates the polynomial coefficients and plots the corresponding function 
on the grid. The program is built in this manner due to limited core require- 
ments within the IBM 1130 computing system. The maximum amount of core 
storage available on this system is 8,000 words of core storage. Since 
GECAP and BESFT each require approximately 8,000 words of core storage, 
it is impossible to execute the program with both routines in core together. 
Execution is achieved through the call LINK system routine which, in effect, 
stores only one routine in core at a time. If the continuous option is used, 
GECAP will call BESFT. When this is done, BESFT is brought into core. 
BESFT executes its function and returns to GECAP. When GECAP returns 
to core only the values for the variables stored in common will be retained. 

Routines: The following is a list of routines called by the mainline 
GECAP. 


1. LABEL - This routine draws and labels the grid and places the 
titles on the graph. 

2. HIST - This routine generates a histogram plot from the input 

data. 


3. ERRO - This routine contains all GECAP error messages and is 
called only when an error is found in the input data. 

4. BESFT - This routine determines the best degree of polynomial 
through a statistical ANOVA test and plots the least squares function on the 
grid. 


5. F - This function calculates the value of the dependent variable 
from the polynomial equation found in the curve fit analysis. 

6. CURF - This routine calculates the coefficients for all degree 
polynomials (up to the 11th degree) which fit the data. 
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7. FINSH - This routine checks to see that all plotting has been 
completed before further execution of the program proceeds. If the plot buffer 
is empty, then the program continues; if not, the program waits for all plotting 
to be completed before continuing. 

NOTE: The routines LABEL, HISTO, and ERRO are placed in a system 
monitor control, LOCAL. This causes these routines to be placed in core 
only when they are called by the mainline GECAP. At all other times these 
routines are stored on the disc. The purpose of this monitor control is to 
prevent the amount of core storage used by GECAP from exceeding 8K words. 


C. Program Restrictions and Limitations 

1. As previously stated, use of a system routine (CALL LINK) and 
a system monitor control (LOCAL) have been made in order to limit the 
amount of core storage required by GECAP. The user should be cautious 
when making program modifications, for these changes could drive core 
requirement over the 1130 computer capacity. 

2. When executing the continuous (C) or histogram (H) options, the 
data must be input in order of increasing x value (x. < x^+1 < x/t-2 < x^+3, etc) . 

3. The statistical test in BESFT does not always guarantee that the 
best degree polynomial has been selected by the program. It is possible that 
the curve selected may be accurate for the given data points but not the best 
"eye balled" curve. This circumstance may be avoided by using more data 
points and by spacing these points evenly with respect to the x axis. 

4. A maximum of 100 data points may be input for any case using 
any option. 

5. The maximum degree equation which will be calculated for any 
case using the continuous option is an 11th degree polynomial. For cases 
where less than 11 data points are input, the maximum degree polynomial 
calculated will be (N-l) where N is the number of data points. 


D. PROGRAM MODIFICATIONS FOR OTHER COMPUTING SYSTEMS 


This program may be modified to run on other systems by making a 
few minor changes on large computers, rather than calling BESFT through 
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a CALL LINK system routine. BESFT should be a subroutine called by 
GECAP and placing a GO TO statement after BESFT sending it to the beginning 
of GECAP. The system monitor, LOCAL, may be deleted. The plotter 
routines must be changed as these routines are unique to each computing sys- 
tem. Also, the input/output unit numbers for all READ and WRITE statements 
must be changed to correspond to the particular installation. 
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APPENDIX A 


DETERMINATION OF THE BEST FIT POLYNOMIAL 


One of the problems associated with plotting out the results of a least- 
squares curve-fit is the determination of the degree of polynomial that best 
fits the empirical data. One possible method is to calculate the least-squares 
coefficients for as many polynomial equations as one desires, and then plot 
out each of the different power curves. Selection of the plotted curves can 
be made against the input data by using the "eye ball’' technique. 7 In many 
cases this is the best method since the analyst has some prior information 
about the shape of the curve, or because higher degree polynomial coefficients 
cannot be determined due to a lack of enough data points. 


However, a routine has been included in GECAP that statistically 
determines the best degree of polynomial to be plotted. The test procedure 
is mathematically defined as follows: 


Given a set of n data points x } , y t ; x 2 , y 2 ; x 3 , y 3 ; + 


x , y ; one 
n J n 


can determine the least squares coefficients b , for an mth degree polynomial, 
assuming m < n-1. An equation of this type is defined in the following form: 


y (x.) = b. + b, x. + b, x? + 
J m v 1 ° 1 i 1 i 


b x 
m i 


m 


The problem is to determine what value of m will give a statistically 
good fit. (i. e. , Does the addition of a higher order term to the polynomial 
equation fit the data significantly better than without using it?) 

A test can be set up that will solve this problem on a probabilistic 
basis. The test uses the ratio of the difference of two independent estimates 
of the error variance against a single estimate of the error variance. Mathe- 
matically, this can be defined as the following F statistic: 

S 2 - S 2 

„ m-1 m 


m 


7. The "eye ball" technique (as named by the author) is simply a method in 
which someone determines the shape of the approximate curve by placing a 
french curve over the given data and estimating what a good fit would look like. 
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where 


X, (y m ( ^- y i )J 

2 = LzJ 

m n - (m + 1) 

is the residual variance after fitting the data with the mth degree polynomial 
and 


.2 , < y m-l (x i> - y i )2 

s2 = 7 T 

m-1 n - ( m) 

is the residual variance after fitting the data with the m-lth degree polynomial. 

The null hypothesis = 0 ( i. e. , the test can be stated as the mth 

degree polynomial is not significantly better than the m-lth degree polynomial 
or the null hypothesis b m = 0) is tested for significance at the 5-percent 

rejection level 
where 


F > F 

m-1 0. 95, v lt v 2 

"l = 1 

v 2 = n-m-1 

to determine if the mth degree curve has a smaller error variance than the 

m-lth polynomial. If the value of F is less than F then it is 

a , v lt v 2 

assumed that the mth degree is satisfactory. This test is also performed 
for the m-2th error variance. After two successive tests show no significant 
in the F values, then it is assumed that the mth degree polynomial is a good 
fit. 8 


It should be pointed out that this test will not always satisfy the user of 
GECAP in all possible cases. However, it works well when the data to be 
fitted is evenly distributed along the x axis and the data has less than two 


8. The corresponding polynomial coefficients are printed out and the mth 
polynomial function is plotted on the graph. 
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inflection points. For data reflecting many (3 or more) inflection points, 
the desired polynomial fit should have at least 30 or more coordinate values 
per curve as input to the GECAP Program. 

More information on this subject can be obtained from the books 
Probability and Statistics for Engineers by I. Miller and J. Freund, page 245, 
and/or Probability and Statistics in Engineering and Management Science by 
W. Hines and D. Montgomery, pages 332 and 356. 
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APPENDIX B 


PROGRAM LISTING 


The following program listing is included for those users wishing to 
modify the existing program to suit a special need. The listing has been well 
commented so that a programmer can isolate any individual logic operation 
within a routine. Information on the IBM 1130/1800 plotter subroutines can 
be obtained from the IBM System Reference Library, Form C26-3755-0. 
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// FOR 

* I OCS ( 1132 PRINTFR.PLOTTER.TyPHWRI TLK.CARJ .KEYBOARD i 
•ONE WORD INTEGERS 
•EXTENDED PRECISION 
C 

C*** GECAP REVISION C* 09/0S/72 
C 


GC APOUUU 
GCAPUOlU 
GCAP0U2O 


£###***#. **#**.-. *####*.*)n***in.*#**#*#*##***#innnnn*.*. *.*.»..-**.-**. 

C GCAPOOAU 

GCAPOUSO 
GCAP0U6U 
GCAP OOVO 
GCAPU08O 
GCAP0090 

s a ****»«gcapoioo 

.CCCCCCCCCLC LGCCCCCC^gLCCCGCAPOI 10 


C*****PR0GRAMMED by 
c 
c 
c 
c 

£**********«**■»»***♦»*»■ 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


B. 

B. 

J. 

R. 


R. BEADLE II 
D. DOLERHIE 
W. OWEN 
A. SCHLAGhECK 




( 9/ /I ) 

( 12/ 71 ) 

19/72) 

(9/71 TO 9/72 ) 

* It * . *#* *. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM GECAP IS A GENERALIZED CURVE — FIT AND PLOT PROGRAM 
FOR USE IN PLOTTING A WIDE VARIETY OF DATA ON A GRID SMALL 
ENOUGH TO FIT ON AN 8 1/2 X 1* SHEET OF PAPER. NUMEROUS 
OPTIONS ARE SUPPLIED Tu THE USER. IN ORDER FUR HIV TO OBTAIN 
A PLOT OF HIS DA (A SET. 


OPTIONS' 


1 “ 


EITHER OR BOTH OF THE AXIS ON THE GRID MAY BE NUMB- 
ERED IN EITHER INTEGER OR REAL NUMBER FORM. 


2 - THE DATA MAY BE PLOTTED IN ANY ONE OF FUUR WAYS 


GCAP012U 
GCAP01 30 
GCAP01 90 
GCAP01 SO 
GCAP016U 
GCAPU1 70 
GCAP 01 80 
GCAP0190 
GCAP0200 
GCAP02 10 
GCAP 02 2 0 
GCAP02 30 


(1) DISCRETELY -THE DATA IS PLOTTED AS I NO I V I DUALGCAP02 AO 

POINTS. WITH THE POINT SYMBOLS DIFFERING FROM GCAP02 SO 
CASE TO CASE. GCAP02 60 

(2) - LINEARLY - THE DATA POINTS ARE CONNECTED bY GCAP0270 

STRAIGHT LINES. GCAP0280 

(31 - CONTINUOUSLY - THE INDIVIDUAL POINTS ARE GCAP0290 

ORIGINALLY PLOTTED AS THE DESCRETE POINTS. BUTGCAP0300 
The POINTS UNDER GO a LEAST-SUUARES CURVE-FITGCAP0310 

AND THE CALCULATED FUNCTION IS PLOTTED. GCAP0320 

(A! - HISTOGRAM - THE DATA POINTS ARE NOT PLOTTED GCAP0330 
INDIVIDUALLY . BUT A HISTOGRAM IS GENERATED 
FROM THE INPUT DATA. 


MULTIPLE CASES MAY BE PLOTTED ON THE SAME GRID WITH- 
OUT INPUTTING ADDITIONAL GRID INFORMATION. THIS 
OPTION DOES NOT ALLOW FOR DIFFERENT TYPE PLOTS ON 
THE SAME GRID. 


IS 


GCAP03AU 
GCAP03 SO 
GCAP0360 
GCAP037U 
GCAP0380 
GCAP03 90 
GCAPOAOU 
GCAPQAlu 
GCAP0A2U 
GCAP0A30 
GCAPOAAO 
GCAPOASO 


-NOTE. ..AN OPTION FOR NO DATA POINT MARKS 
AVAILABLE FOR ANY TYPE OF PLOT. 

-NOTE IF DATA POINT MARKS ARE USED. AN OPTION MAY 

BE USED TO LABEL EACH CURVE INDIVIDUALLY. 

-NOTE IF THE USER DESIRES TO SHOW THAT HIS LABELED GCAP0A60 

VALUES ARE RAISED TO A POWER OF 10. AN OPTION EXISTS GCAP0A70 
WHICH LABELS THE LAST 'TIC' MARK ON THE AXIS IN SUCH GCAP0980 
A MANNER TO INDICATE HOW THE VALUES SHOULD BE READ. GCAP0490 

GCAP0300 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCGCAPOS 10 


INTEGER TYPE ( 3 ) . DESC 
INTEGER CONT 
INTEGER P.R 
INTEGER HISTO.BLNK 
INTEGER O.ERR 


GCAPOS 20 
GCAP05 30 
GCAP0590 
GCAPOS 50 
GCAP0S60 
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nn n non non nonn 


INTEGER TEN ( 2 ) 

DIMENSION XNAMEI5) 

COMMON T I TLX ( 7 ) , T I TLY ( 7 ) »T I TLL ( 7 ) 
COMMON XVALUI 101 ) .YVALU! 101 ) 
COMMON <1 .XMAX. YMAX.X INC tXSTRT 
COMMON YSTRT .XSCLE .YSCLE 
COMMON NCASE.ICASE 
COMMON \OCAT 


GCAP06 / J 
GC AP05 80 
GCAP05 9u 
GCAP0600 
GCAPOa 10 
oC AP06 2 o 
GL AP Ofa J O 
GLAP06AU 


DATACONT/ 'C 1 / 
DATADESC/'O' / 

OAT AlINE/ ' L' / 
DATAHISTO/ 'H' / 
DATAXCHV/'END'/ 

DAT ABLNK/ ' •/ 

DATABLANK/' '/ 

OATAACK/’OF '/ 

ICODE ■ 0 
0 = 0 
K = 0 
R = 2 
P * 3 
TYPE ( 3 ) =1 
IDUMY = LINE 


GCAP 06 D U 
GLAP066O 
GC AP Oo 7 o 
UCAP0680 
GLAP069 u 
GCAP J /OO 
OCAP071O 
GCAPO /20 
GCAPG7 30 
GCAP07A0 
GCAPO 760 
GCAPO /60 
GCAPO 7 70 
GCAP07 80 
GCAPO 7 90 


10 

20 

30 

AO 

60 

60 


READ IN THE INFORMATION ABOUT 
THE CURVES. AND THE EXPONENTS 


THt 

FOR 


NUMBER OF 
THE X AND 


CASES . 
Y AXES 


THt LABtLS FOR 


READ(R.20)TITLZ.ICASE.A.NCASt. IXNAME(O) .0= 1 . 5 ) . TtN t 1 ) . TEN ( 2 ) 
FORMAT ( AA.1X.I1. IX. A3. 11. 3X.5A6.3X.I2.3X. 12) 

IF (NCASE-1 1A0.30.30 
IF ( A-ACKI AO. 60, AO 
ERR = 1 

CALL ERROI TYPE, ERR, K.YINC) 

CONTINUE 

IF ( I CASE-1 170.70.60 
IF INCASE- I CASE 18 20,3 70,3 70 


*** THIS IS THE CONSOLE MESSAGE FOR THt USER 


GCAP0800 
GCAPOd 10 
GCAP Od 2 0 
GCAP 08 30 
GCAPOd AO 
GCAP085O 
GCAP 0860 
GCAPOd 7U 
GCAPoaao 
GCAPOB 90 
GCAP090U 
GCAP091O 
GCAP0920 
GCAP0930 
GCAP09A0 


GCAPQ960 

70 WRITE! 1.801 GCAP0960 

80 FORMAT! ' **** POSITION PLOTTER PEN APPROXIMATELY 3 INCHES FROM THtGCAP0970 
1 ***#.,/,< **** RIGHT EDGE OF The PLOTTER. . . . PRESS START . GCAP0980 
2, , *##*',//) GCAP0990 

PAUSE GCAPloOU 

GCAP10 10 

*** READ IN THE INFORMATION ABOUT THE GRID GCAP1020 

GCAP1U30 


READ! R.90 > X I NC » Y INC »XMAX .YMAX.XSTRT • YSTRT .TYPE ! 1 > .TYPE ( 2 I .TYPE I 3 ) .GCAP 10 AO 
1N0DAT GCAP 1050 


90 FORMAT (6F10.0.AA1) 

'** PRINT OUT THE INPUT DATA 

WRITE(P.IOO) 

100 FORMAT ( lHl ) 

WRITEIP.llOIXINC.YINC .XMAX.YMAX.XSTRT .YSTRT .NCASE 
110 FORMAT!/, • X-AXIS INCREMENT = 1 »F 10. 2 »5X * ' Y-AXI S I N< 
1,/.' X-AXIS LIMIT = 1 »F10. 2. 5X. 'Y-AXIS LIMIT = 

2AXIS ORIGIN = ' .F10. 2 ,5X. ' Y-AXIS ORIGIN = 'F10 

3 = '.110) 


GCAP 1060 
GCAP 1070 
GCAP 1080 
GCAP 1090 
GCAP1100 
GCAP 1 1 10 
GCAP 1120 

IREMENT =• .F10.2GCAP1130 
•F10. 2,/,' X— GC AP 1 1 AO 
,2»/»' NCASES GCAP 1150 

GCAP 1160 
GCAP 1170 
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U V V U V U V 


C*** CHECK THE VALUES FOR LABELING T HE AXES TO SEE IF T ME STARTING 
C*** POINTS ARE LESS THAN T Ht END POINTS 
C 

IF ( XMAX-XSTRT H2D.120.130 
120 ERR*2 

CALL ERRO(TYPE.CRR.K.YINC) 

GO TO 150 

130 IF! YMAX-YSTRT > 140,140. 160 
140 ERR-3 

CALL ERRO(TYPE.ERR*K«YINC) 

150 STOP 
160 ERR = 4 

CALL ERRO(TYPE*ERR.K*YINC) 

C 

C#** CONVERT ALPHA-NUMERIC INPUT DATA TO NUMERIC VALUES 
C 

DO 260 K» 1 . 3 

WRITE ( P* 170IK *TYPE (K ) 

170 FORMAT!' TYPE!'. 12.') = ’.6X.AA) 

IF ( TYPE (K l-CONT ) 190 . ISO . 190 
180 TYPE ( K ) = 1 
GO TO 260 

190 IF ( TYPEIKl-DESC 1210*200.210 
200 TYPE ( < ) = 2 
GO TO 260 

210 IF (TYPE(K)-HISTO)230. 220.230 
220 T YPE I K ) =4 
GO TO 260 

230 IF (TYPE(K)-LINE)250.240.250 
240 TYPE ( < ) =3 
GO TO 260 
250 ERR-5 

CALL ERRO(TYPE.ERR.K.YINC) 

260 CONTINUE 
C 

C*** READ IN TITLES FOR THE X AND Y AXIS AND ALSO THE GRAPH TITLE 
C 

READ! R .270 )(TITLX(I).I=1.7),(TITLY(J).J=1.7).(TITLE(K),K=1.7) 
270 FORMAT (6A4.A2.6A4.A2. 7 A 4) 

280 CONTINUE 


GCA ’ . 1 So 
UCAPilSv, 
GCAP 12GO 
GCAP12 10 
GCAP 1220 
GCAP 12 30 
GCAP 1240 
GCAP 1250 
GCAP126U 
GCAP 12 70 
GCAP 12 60 
GCAP 1 2 90 
GCAP 130U 
GCAP 13 lu 
GCAP 13 20 
GCAP 1330 
GCAP 1340 
GCAP1350 
GCAP 1360 
GCAP 13 70 
GCAP 1380 
GCAP 1 3 9 0 
GCAP 1 4u0 
GCAP 14 1 U 
GCAP142U 
GCAP 1430 
GCAP 1 440 
GCAP 1450 
GCAP 1460 
GCAP 14 70 
GCAP 14 8 0 
GCAP 1490 
GCAP 1300 
GCAP 15 10 
GCAP1520 
GCAP 1530 
GCAP 1540 
GCAP1350 
GCAP 1560 
GCAP 15 70 


THIS SECTION COMPUTES THE VALUES FOR THE GRID FROM THE SPECIFIED 

USER INFORMATION. 


INCX - X I NC*1 0* 

INCY = Y I NC*1 0 • 

1X1 - INCX/10 
I Y I - INCY/10 

IF! ( IYI*1 0 1- INCY 1300.290. 300 
290 IF! ( 1X1*10 ) -INCX ) 310.320 .310 
300 I Y I » IYI+1 
GO TO 290 
310 IXI - IXI + 1 
320 NXTCS * (XMAX-XSTRT l/XINC 

THIS IF TEST LIMITS THE NUMBER OF 

!F(NXTC5-20)340. 340.330 
330 ERR«6 

CALL ERROITYPE.ERR.K.YINCI 


GCAP 1590 
GCAP 1600 
GCAP 16 1U 
GCAP 1 6 2 0 
GCAP163U 
GCAP 1640 
GCAP 16 50 
GCAP 16 6U 
GCAP 1670 
GCAP 1680 
GCAP 16 90 
GCAP 1 7 00 
GCAP17 10 
GCAP 1720 
GCAP 17 30 

'TIC' MARKS FOR THE X-AXIS GCAP1740 

GCAP17 50 
GCAP 1760 
GCAP 1770 
GCAP1780 


\ 
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GO TO 280 

340 NYTCS * I YMAX-YSTRT l/YINC 
C 

C THIS IF TEST LIMITS THE NUMBER OF 'TIC' MARKS FOR THE Y-AXIS 

C 

I F ( NYTCS - 20 1360. 360 .3 50 
350 ERR-7 

CALL ERROITYPE. ERR. K.YINC) 

GO TO 280 
360 CONTINUE 

CALL LABEL I T YPE .TEN. IX I . I YI . I NCX . INCY .NXTCS .N YTCS . YI NC > 

370 CALL SCALE ( XSCLE .YSCLE .XSTRT .YSTKT ) 

C 

c*** READ IN DATA POINTS AND PLOT THEM 
C 

K I ■ 1 

CALL EPLOT ( 1. XSTRT .YSTRT ) 

3 80 RE AD l R.3 90) XV ALU (K I I .YVALUIKI ) .XCH 
390 FORMATI2F10.3.A?) 

IF (XCH-XCHV 1400 .490.400 
400 IF (YVALUIKI ) -YSTRT ) 440 .4 10 .41 0 
410 IF (YVALUIKI ) -YMAX I 420 . 420 . 440 
420 IFIXVALUIKI ) -XSTRT 1440.430.430 
430 IFIXVALUIKI ) -XMAX ) 450 .450 . 440 
440 ERR-9 

CALL ERROI TYPE. ERR. K.YINC) 

450 CONTINUE 

460 IF I KI — 101 1470.480.480 
470 K I ■ K I + 1 
GO TO 380 
C 

C**# A LIMIT OF 100 DATA POINTS IS ESTABLISHED FOR THIS PROGRAM 
C 

480 ERR*8 

CALL ERROITYPE. ERR. K.YINC) 

GO TO 380 
490 NCK-TYPEI3) 

GO TO (500.510.510.830) .NCK 
500 JFK - 1 
GO TO 520 
510 JFK • 2 

520 GO TO(540.530)»JFK 
530 ICODE - 2 
540 JG-KI-1 



IF I TYPE I 3 >-3)550.560.550 
550 IFINODAT-BLNK1730.560.730 
560 DO 630 LBJ-l.JG 

IF I TYPE 1 3 1-3)600.570.600 
570 IFILBJ-l 1600.580.590 
580 CALL EPLOTI l.XVALUILBJ) .YVALUILBJ) ) 

IF I NODAT-BLNK 16 30.610.630 
590 CALL EPLOTI2.XVALUILBJ) .YVALUILBJ) ) 

IF I NODAT-BLNK 1630 .610 >630 
600 CALL EPLOTI 1 .XVALU I LB J ) . YVALU I LB J ) ) 

CALL E PLOT I 2. XVALU I LB J) .YVALUILBJ) ) 

610 1 1 » I CASE-1 
620 CALL POINT! II ) 

630 CONTINUE 
C 

C*** LABEL THE CURVE 


GCAP1790 
GC AP 18 00 
GCAP1810 
GCAP182U 
GCAP 1830 
GCAP 1840 
GCAP 1850 
GCAP I860 
GCAP1870 
GCAP 1880 
GCAP 18 9 J 
GCAP 1900 
GCAP19 10 
GCAP 19 20 
GCAP1930 
GCAP 1940 
GCAP1950 
GCAP1960 
GCAP1970 
GCAP1980 
GCAP1990 
GCAP2000 
GCAP20 10 
GCAP2020 
GCAP20 30 
GCAP2040 
GCAP2050 
GCAP2060 
GCAP2070 
GCAP2080 
GCAP2090 
GCAP2100 
GCAP2110 
GCAP2120 
GCAP2130 
GCAP2140 
GCAP2150 
GCAP2160 
GCAP2170 
GCAP2180 
GCAP21 90 
GCAP2200 
GCAP22 10 
GCAP2220 
■GCAP2230 
GCAP2240 
6CAP2250 
GCAP2260 
GCAP2270 
GCAP22 80 
GCAP2290 
GCAP2300 
GCAP2310 
GCAP2320 
GCAP2330 
GCAP2340 
GCAP2350 
GCAP2360 
GCAP2370 
GCAP2380 
GCAP2390 
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DO 640 0=1 ,5 

IF (XNAME(O)-ULANX) 6 50.640.6 50 
640 CONTINUE 
GO TO 730 

650 X= ( XMAX )-( XMAX-XSTRT ) / (4.5 ) 

660 Y=YSTRT-(.5)/ ( YSCLF ) -(.2)*(II)/(YSCLE) 

670 CALL EPLOT (l.X.Y) 

CALL EPLOT ( 2 * X * Y ) 

680 CALL POINT (II) 

690 X= ( XMAX )-< XMAX-XSTRT ) /( 5. ) 

Y*Y-( . 05 ) / ( YSCLE ) 

700 CALL ECHAR ( X * Y . . 1 . . 1 .0 . 0 ) 

710 WRITE (7.720) (XNAME(O) .0=1.51 
720 FORMAT ( 5A6 ) 

C 

730 WRITEIP.740) TITLX.TITLY 

740 FORMAT (///. 14X . ' ***INPUT DAT A*»* ' . / / . 1 X . 6 A4 • A2 . 4X .6A4 . A2 » / / 
WRITE ( P.750 ) T ITLZ. I CASE. A. NCASE 
750 FORMAT(20X.A4,I2.1X.A2.I2) 

760 XI = XI - 1 
DO 780 I = 1 . < I 

WR I TE ( P » 770 ) XVALU ( I ) . YVALU ( I ) 

770 FORMATI7X.F10.4.14X.7F10.4) 

780 CONTINUE 

XI « XI + 1 

CALL EPLOT ( 1 .XSTRT .YSTRT ) 

790 CONTINUE 

I F ( I CODE-2 1800.810.810 
800 CALL FINSH 

CALL LI NX ( BESFT ) 

810 IF(NCASE-ICASE>820.820,10 
C 

C*** SET THE PEN FOR THE NEXT PLOT 
C 

820 CALL EPLOT ( 1 .XMAX. YSTRT ) 

CALL SCALE! 1.0. 1.0. 0.0. 0.0) 

CALL EPLOT ( 1 .4. 0.0.0 ) 

CALL EXIT 
830 CALL HIST 
ICODE-2 
GO TO 730 
END 

// DUP 

♦DELETE GECAP 

♦STORE WS UA GECAP 


GCAP240 J 
GCAP2410 
GCAP242 J 
GC AP 24 3 0 
GC AP 244 0 
GC AP 24 5 0 
GCAP246 -j 
GC AP 2<* 79 
GCAP2480 
GCAP2490 
GC AP 2 5 00 
GCAP25 10 
GC.AP252 J 
GC AP 2 5 3 0 
GCAP2540 
•GCAP2550 
GCAP2560 
GCAP2570 
GCAP2 6 80 
GCAP2590 
GCAP2600 
G CAP 2 6 10 
GCAP2620 
GCAP2630 
GCAP2640 
GCAP2650 
GCAP2660 
GC AP 2 6 70 
GC AP 2 6 8 0 
GCAP269Q 
GC AP 2 7 0 0 
GCAP27 10 
GCAP27 20 
GCAP2730 
GC AP 2 740 
GCAP2750 
GCAP2760 
GCAP2770 
GCAP2780 
GC AP 2 7 9U 
GCAP2800 
GCAP28 10 
GCAP28 20 
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// FOR 

♦ONE WORD INTEGERS 
•EXTENDED PRECISION 

C LABLOOOO 

C*«* THIS SUBROUTINE DRAWS Tt-.E GRID. LABELS THE GRID. AND LABL0013 

C*** PLACES THE TITLES IN POSITION. LABL002O 

C LABLOO 30 

SUBROUTINE LABEL (TYPE .TEN. I XI . I YI . INCX. INCY .NXTCS.NYTCS.YINC ) LABLOO AO 

C LABL00SO 


C*** 

c 


c 

c*** 

c*** 

c*** 

c*** 


GECAP REVISION C* 09/05/72 

INTEGER TYPE ( 3 ) . TEN ( 2 ) 

COMMON TITLXI 7) .T I TLY ( 7 ) .T I TLE < 7 ) 

COMMON XVALUI 101 ) .YVALU( 101 ) 

COMMON KI.XMAX.YMAX.XINC. XSTRT 
COMMON YSTRT.XSCLE.YSCLE 
COMMON NCASE. ICASE 
COMMON NODAT 

THE LENGTH OF THE X AND Y AXES ARE GOVERNED BY THE 
MENTS . FOR LONGER AXIS. CHANGE THE NUMBERS IN THE 
-MENTS SO THAT THE ACTUAL LENGTH OF THE AXIS UN I 

in the dividend of the statement. 


LA9L0060 
LABL0070 
LABLOOH J 
LABL0090 
LABL0100 
LABL01 10 
LABL012U 
LABL013U 
LAHL01 AO 
LABL0150 
NEXT TWO STATE-LABL0160 
FOLLOWING ST ATELABL01 70 
NCHES ) WILL BE LABLOi 80 
LABL0190 


C 


C 

c 

c*** 

c 


c 


XSCLE • 8 « 5/ ( XMAX - XSTRT ) 
YSCLE • 6.0/< YMAX-YSTRT) 
YINT - “.2 /YSCLE 
XINT » -.32/XSCLE 


DRAW THE GRID 

CALL EPLOT ( 1.0. .0. ) 

CALL EPLOT ( 2 .0. • 0. ) 

CALL SCALEtXSCLE. YSCLE. 0. .0. ) 
CALL EGRIDt 1.0..0..YINC. NYTCS) 
CALL EPLOT ( 1 .0. .0. ) 

CALL EPLOT ( 2 .0. .0. ) 

CALL EGR ID ( 0 .0. . 0 . .XI NC * NXTCS ) 


C*»* LABEL THE GRID 
C 

NXTCS ■ NXTCS ♦ 1 
NYTCS • NYTCS + 1 
X • XINT 
XNUM « XSTRT 
INIT ■ XSTRT 
NUM ■ XSTRT«10. 

IF ITY PE (3) “A >20.10.20 
10 X-(XINTI+(XlNC>/(2.) 

NXTCS-NXTCS-1 
20 DO 80 I “1 .NXTCS 

CALL ECHARIX. YINT. .1.. 1.0.0) 
IFlTYPE(l)-2>60.30.60 
30 WRITEI7.A0) INIT 
AO FORMAT (14 1 
50 FORMAT I I 5 > 

INIT-INIT+IXI 
GO TO 70 

60 WRITEI7.520) XNUM 
70 CONTINUE 


LABLC200 
LABL02 10 
LABL0220 
LABLC2 30 
LABL02AG 
LABL0250 
LABL026U 
LABL0270 
LABL02 80 
LABL02 90 
LABL0300 
LABL03 10 
LABL0320 
LABL0330 
LABL03A0 
LABL0350 
LABL0360 
LABL0370 
LABL0380 
LABL0390 
LABLOAOO 
LABLOA 1 0 
LABL0A20 
LABL0A30 
LABLOAAO 
LABL0A50 
LABL0A60 
LABLOA 70 
LABL0A80 
LABL0A90 
LABL0500 
LABL0510 
LABL0520 
LABL05 30 
LABL05A0 
LABL0550 
LABL0560 
LABL0570 
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n r> n 


X=X*XINC 

NUM-NUM+INCX 

XNUM * FLOATtNUM ) /10. 

BO CONTINUE 

*** PLACE EXPONENT AT THE END OF THE X AXIS 

IF ( TEN! 1 ) 190. 270*90 
90 IFITYPEI 1 1-2)1,70.100.170 
100 X*XMAX+.2/XSCLE-XSTRT 
Y«YINT 

IF(TYPE(3>-4) 120.110.120 
110 X«X-XINC/2. 

120 CALL ECHAR (X.Y..1. .1.0.0) 

130 WRITE (7.300) 

X-X+.4/XSCLE 
I F ( TEN ( 1 ) 1160.160.140 
140 IF(TENll)— 10)150.160*160 
150 X 3 X— • 1 /XSCLE 
160 Y=Y+.05/YSCLE 

CALL ECHAR ( X .Y . . 1 . . 1 .0. 0 ) 

WRITE (7*310)TEN(1) 

GO TO 270 

170 X=XMAX+.4/X5CLE-XSTRT 
Y-YINT 

IF (TYPE (3 1-4)190.180.190 
180 X-X-XINC/2. 

190 CALL ECHAR ( X »Y . . 1 . . 1 .0.0 ) 

200 WRITE (7.300) 

210 X-X+.4/XSCLE 

I F ( TEN( 1 ) 1240.240.220 
220 IF(TEN(1)-10)230.240.240 
230 X=X”« 1 /XSCLE 
240 Y*Y+.05/YSCLE 
250 CALL ECHAR (X.Y.. It. 1.0.0) 

260 WRITE (7*310) TEN ( 1 ) 

270 CONTINUE 
C 

C*** LABEL THE Y AXIS AND PLACE EXPONENT AT THE END OF THE Y AXIS 
C 

XINT — .6/XSCLE 
IF (TYPE( 2 1-2)280.290.280 
280 XI NT "-.6 5 /XSCLE 
290 Y-0.0 

INIT * YSTRT 
NUM - YSTRT 
XNUM" YSTRT 
DO 540 I "1 .NYTCS 
300 FORMAT ( 'X 10' ) 

310 FORMAT! 12) 

320 IF(NYTCS-I 1490.330.490 

330 I F ( TEN < 2 J 1340.490.460 

340 FF"0.0 

350 X*— .7/XSCLE+FF 

360 CALL ECHAR’ (X.Y..1. .1.0.0) 

370 WR I TE ( 7 . 300 ) 

380 X — .3/XSCLE 

390 ZZ"Y+.05/YSCLE 

400 CALL ECHAR(X.ZZ».l.. 1.0.01 

410 WRITE(7*310)TEN(2) 

420 X I NT* XI NT-.7/XSCLE+FF 


LABL0580 
LABL05 90 
LABL0600 
LABL0610 
LABL06 2 0 
LABL0630 
LABL0640 
LABL0650 
LABL0660 
LABL06 70 
LABL0680 
LABL06 90 
LABL070 o 
LABL071U 
LABL0720 
LABL0730 
LABL0740 
LABL075 0 
LABL0760 
LABL0770 
LABL0760 
LA8L0790 
LABL0800 
LABLOB 10 
LABL0820 
L ABL 08 30 
LABL0840 
LABL0850 
LABL0860 
LABLO870 
LABL0880 
LABL0890 
LABL0900 
LABL0910 
LABL0920 
LABL0930 
LA6L0940 
LABL0950 
LABL0960 
LABL0970 
LABL09 80 
LABL0990 
LA8L1000 
LABL1010 
LA8L1020 
LABL1030 
LABL1040 
LABL1050 
LABL1060 
LABL1070 
LABL1080 
LABL1090 
LABL1100 
LABL11 10 
LABL11 20 
LABL1130 
LABL1140 
LABL1150 
LABL1160 
LABL1170 
LABL1180 




430 IF ( TYPE ( 2 ) -2 1440 ,490.440 
440 XINT«XINT-.05/XSCLE 
450 GO TO 490 

460 IFIT£N(2)-iO)470.340,340 
470 FF=.1/XSCLE 
480 GO TO 350 
490 CONTINUE 

CALL ECHARIXINT.Y..1. .1.0.0 
IF ( TYPE! 2 1-2 ) 510*500.510 
500 WRITEI7.50) INIT 
INIT = INIT +IYI 
GO TO 530 

510 W R I T E ( 7 * 5 2 0 ) X MUM 
520 FORMAT ( F6* 2 ) 

530 CONTINUE 

Y = Y + YINC 
NUM * NUM + INCY 
XNUM ■= XNUM + YINC 

540 CONTINUE 

*** PLACE THE TITLES ON THE GRAPH 

X => 2.4/XSCLE 

Y = -l./YSCLE 

CALL ECHARIX.Y,. 15. .15*0.0) 
WRITEI7.550) ( T I TLX ( L ) . L = 1 . 7 ) 
550 FORMAT ( 7A4 ) 

ANG » 3.14157/2.0 
X « -.68/XSCLE 
I F ( TYPE (2 >-2 ) 5 60 .5 70 .560 
560 X — .73/XSCLE 
570 Y-1.0/YSCLE 

CALL ECHARIX.Y. .15. .15. ANG) 
WRITE! 7.550) (TITLY(M) ,M= 1,7) 

X « 2. 5/XSCLE 

Y = 6. 5/YSCLE 

CALL ECHARIX.Y. .2. .2, 0.0) 
WRITE (7. 550) (TITLE(N) .N»1.7) 
CALL EPLOT ( 1 .0. ,0. ) 

CALL EPLOT I 2 • 0. , 0 . ) 


RETURN 


END 


// DUP 


•DELETE 

LABEL 

•STORE WS 

UA LABEL 


LAHLI190 
LABL120U 
LABL12 10 
LABL1220 
LABL123J 
LABL1240 
LA BL 12 50 
LABL1260 
LABL127J 
LABL12B0 
LABL129U 
LABL1300 
LABL1310 
LABL1320 
LABL1330 
LABL 1340 
LABL1350 
LABL 1360 
LABL 1 3 70 
LABL 1 3 8 0 
LABL 1390 
LABL1400 
LABL 1410 
LABL 1420 
LABL 1430 
LABL 1 44 U 
LABL1450 
LABL 1460 
LABL1470 
LA8L1480 
LABL 1 490 
LABL 1 500 
LABL15 10 
LABL1520 
LABL1530 
LABL 1 540 
LABL 1 5 50 
LABL 1 560 
LABL 1 5 70 
LABL 1580 
LABL 1590 
LABL 1600 
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// FOR 
•ONE WORD 
•EXTENDED 
C 


INTEGERS 

PRECISION 


C*** 

c*** 

c 

c 

c*** 

c 


THIS SUBROUTINE CONTAINS ALL 
CALLED ONLY BY GECAP 


ERROR MESSAGES FOR GECAP AND IS 


SUBROUTINE ERRO ! T YPE . ERR .X . Y I NC ) 

GECAP REVISION C. 09/05/72 

INTEGER TYPE! 3) .ERR 

COMMON TITLX(7)»TITLY(7I »T I TLE I 7 ) 

COMMON XVALU! 101 ) .YVALU! 101 ) 

COMMON XI .XMAX.YMAX.XINC.XSTRT 
COMMON YSTRT.XSCLE.YSCLE 
COMMON NCASE.ICASE 
COMMON NODAT 

GO TO ( 10.30.50.70.210.250.270.280.300) .ERR 


10 WRITE (1.20) 

20 FORMAT! /•' ****ERROR**«* CASE CARD INPUT 
l.'*NUMBER OF CASES HAS BEEN SET EQUAL TO 
NCA5E* 1 
GO TO 320 
30 WRITE! 1.40) 

40 FORMAT! /. ' ****ERROR****END POINT OF X 
10INT****' ./»• ****EXECUT I ON 

GO TO 320 
50 WRITE! 1.60) 

60 FORMAT! /« ' ****ERROR****END POINT OF Y 
10INT****' ./. ' ****EXECUT ION 

GO TO 320 


VARIABLE IS INCORRECT ' 
1 + EXECUTION RESUMED*' 


AXIS IS LESS THAN STARTING 
DI SCONT I NUED**** ' ) 


AXIS IS LESS THAN STARTING 
DISCONTINUED****' ) 


C 

c#** 

c*** 

c*** 

c 

70 


CHECK THE NUMBER OF DECIMAL PLACES IN 
THE AXES AND WRITE OUT AN ERROR MESSAGE 
DECIMAL PLACES ARE USED 

NCXI-IFIXIXINC) 

CNX-IXINC-NCXI )*<100. ) 

NCX-IFIX(CNX) 


THE VALUES USED FOR LABEL I 
IF MORE THAN TWO 


ERROOOOO 
ERROOO 1 0 
ERR00020 
ERR00030 
ERR00040 
ERR0005U 
ERR00060 
ERR0007 0 
ERR00080 
ERR00090 
ERR00100 
ERR001 10 
ERR00120 
ERR001 30 
ERR001 40 
ERR001 50 
ERR00160 
./ERR001 70 
> ERR001 80 
ERR001 90 
ERR002 00 
ERR002 10 
PERR002 20 
ERR002 30 
ERR00240 
ERR00250 
PERR002 60 
ERR002 7u 
ERR002 8 0 
ERR00290 
NGERR00300 
ERR003 1 0 
ERR00320 
ERR003 30 
ERR00340 
ERR003 50 
ERR00360 


IF(CNX-NCX)80.90.80 
80 WRITE! 1 .190) 

PAUSE 

90 NCYI-IFIX! YINC) 

CNY-(YINC-NCYI)*(100. ) 
NCY* IFI X ( CNY ) 
IFICNY-NCYI100.110.100 
100 WRITE! 1.190) 

PAUSE 

110 NXX* IF I X I XSTRT ) 

XXX-(XSTRT-NXX)*(100. ) 
I XX* IF I X ( XXX ) 

IF (XXX-I XX) 120.130.120 
120 WRITE! 1.190) 

PAUSE 

130 NYY-IFIX! YSTRT) 

YYY-(YSTRT-NYY)*(100. ) 
IYY-IFIX! YYY) 

IF! YYY-IYY) 140.150.140 
140 WRITEI1.190) 

PAUSE 


ERR003 70 
ERR003 8 0 
ERR00390 
ERR00400 
ERR004 10 
ERRQ042 0 
ERR00430 
ERR00440 
ERR004 50 
ERR004 60 
ERR00470 
ERR00480 
ERR00490 
ERR005 00 
ERR005 10 
ERR00520 
ERR00530 
ERR00540 
ERR00550 
ERR00560 
ERR005 70 
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150 NXM-IFIX(XMAX) 

XXM» ( XMAX— NXM ) * ( 100. ) 

I XM* I F I X ( XXM ) 
IFIXXM-IXMI160.170.160 
160 WRITE! 1.190) 

PAUSE 

170 NYM-IFIX(YMAX) 

YYM« ( YMAX— N YM ) * ( 100 • ) 

I YM = IF I X 1 YYM I 
IF< YYM-IYM >180.200.180 
180 WRITE! 1.190) 

PAUSE 


ERR005 80 
ERR00590 
ERR00600 
ERR006 1 0 
ERR006 20 
ERR006 30 
ERR00640 
ERRO065O 
ERR00660 
ERR00670 
ERR006R0 
ERR006 90 


190 FORMAT(2X.'****ERROR****THE VALUES SPECIFIED FOR LABELING THE AXESERR00700 

1' ./.'REQUIRE MORE PRECISION THAN THAT ALLOWABLE IN GECAP ERR00710 

2 ' »/»'•«. CHECK GECAP USER MANUAL FOR DET A I LS . . . ' . / . ' TH I S ERROR MAY ERR0072U 
3BE INSIGNIFICANT IN SPECIAL CASES. ..IF EXECUT I ON ' . / . ' I S STILL DESI ERRC07 30 
4RED PRESS START... IF NOT. PRESS STOP') ERR00740 

200 GO TO 320 ERR00750 

210 WRITE! 1.220IK ERR00760 

220 FORMAT ( 2X . ' ***#ERR0R**** INCORRECT PLOT OPTION WAS USED FOR'./.' ERR00770 

1 VAR I ABLE TYPE t ' . 1 1 . ' ) .****CHECK GECAP USER MANUAL FOR DETAILS ERR007S3 

2 ' ./ . ' ****TYPE IN THE DESIRED OPTION FROM THE CONSOLE TYPEWRITER ERR00790 
3' ./.' IC.D.H. OR L). ****PRESS (EOF) BUTTON TO CONTINUE EXECUT ION*ERRO080u 


4#*## i ) 

230 READ ( 6 .240 I TYPE ( X ) 

240 FORMAT ( A 1 ) 

K-K-l 
GO TO 320 
250 WRITE! 1.260) 

260 FORMAT!/.' ***ERROR*** 


ERR008 10 
ERR0082U 
ERR0083U 
ERR00840 
ERR00850 
ERR00860 

NUMBER OF TIC MARKS EXCEEDS LEGIBILITY LIERR00B7Q 


1MI T *** ' . / . ' STANDARD FIX-UP TAKEN. EXECUTION CONTINUING') ERR00880 

NXINC - XMAX/10. ERR008 90 

XINC ■ NXINC ERR00900 

GO TO 320 ERR009 10 

270 WRITE! 1.260) ERR00920 

NYINC * YMAX/ 10 . ERR00930 

YINC « NYINC ERR00940 

GO TO 320 ERR00950 

280 WRITE! 1.290) ERR00960 

290 FORMAT!/.' **** ERROR. . MORE THAN 100 DATA POINTS ON INPUT **** ' .ERR009 70 
l/,t **** EXECUTION CONTINUING WITHOUT REMAINING POINTS**') ERR00980 

GO TO 320 ERR00990 

300 WRITE! 1.310) ERROIOOO 

310 FORMAT ( 2X • ' ****ERROR****A DATA POINT WAS FOUND TO EXCEED THE L 1 M I TERR010 10 
IS OF THE AXES...****'./.' ****EXECUT I ON DISCONTINUED****') ERR01020 

GO TO 330 ERR01030 

320 RETURN ERR01040 

330 STOP ERR01050 

END ERR01060 

// DUP 

•DELETE ERRO 

•STORE WS UA ERRO 
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n n n 


// FOR 

*ONE WORD INTEGERS 
•EXTENDED PRECISION 


*** THIS SUBROUTINE GENERATES A HISTOGRAM PLOT FROM THE INPUT DATA 


SUBROUTINE HIST 
C 

C*** GECAP REVISION C. 09/05/72 

C 

COMMON TITLX(7),TITLY<7).TITlE< 7) 
■COMMON XVALUI 101 ) .YVALUl 101 ) 

COMMON K I .XMAX.YMAX.XINC tXSTRT 
COMMON YSTRT .XSCLE * YSCLE 
COMMON NCASE « I CASE 
10 NPTS-KI-1 
20 Z=XVALU(2)-XVALU(1) 

30 IFIXVALUI 1 I-XSTRT 1230.60. AO 
AO CALL EPLOT < 1 .XVALU ( 1 > .YVALU ( 1 ) > 

50 CALL EPLOT ( 2 .XV ALU 111 . YS I R T > 

60 DO 100 I * 1 * NPTS 

70 CALL EPLOT ( 1 .XVALU II) .YVALU ( I > > 

80 PLUS* X VALU ( I ) +2 
90 CALL EPLOT ( 2 .PLUS . YVALU ( I ) I 
100 CONTINUE 

110 CALL EPLOT < 1 . XS TRT . YSTRT I 
120 NPTS-NPTS-1 
130 DO 200 1*1 *NPTS 

1 AO IF (YVALUtI )— YVALU (I+ll 1180*150*150 
.150 CALL EPLOT ( 1 .XVALU ( I +1 ) « YVALU ( I ) ) 
16.0 CALL EPLOT ( 2 .XVALU ( I +1 > • YSTRT ) 

170 GO TO 200 

180 CALL EPLOT ( 1 *XVALU I I +1 ) * YVALU ( I + 1 ) ) 
190 CALL EPLOT ( 2 tXVALU ( I +1 > .YSTRT ) 

200 CONTINUE 

210 CALL EPLOT ( 1 .PLUS * YVALU ( I) ) 

220 CALL EPLOT ( 2 .PLUS . YSTRT ) 

230 RETURN 
END 

// DUP 

•DELETE HIST 

•STORE WS UA HIST 


histoooo 

HI ST 00 10 
H I ST 0020 
HIST 00 30 
HISTOOAO 
H I ST 00 5 u 
H I S TOO 6u 
HI5T0070 
HIST0080 
HIST0090 
HIST0100 
HIST0110 
HIST0120 
HIST0130 
HIST01A0 
HIST0150 
HIST0160 
HIST0170 
HIST 0180 
HIST 0190 
HIST 0200 
rt I ST02 10 
HIST0220 
HI ST0230 
HIST02A0 
HIST0250 
HIST0260 
HIST0270 
HIST0280 
HIST0290 
HIST0300 
HIST0310 
HIST0320 
H I ST0330 
HIST 03 AO 
HIST0350 
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// FOR 

*IOCS< 1132 PR INTER. KEYBOARD .TYPEWR I TER .CARD. PLOT TER ) 

•EXTENDED PRECISION 
•ONE WORD INTEGERS 

C BSTFOOOU 


C*** 

c 

c 

c 

c 

c 


GECAP REVISION C. 09/05/72 


PROGRAM 8ESFT WAS DESIGNED TO SUPPLY A BEST-FIT LEAST SUUARES 
FUNCTION TO A GIVEN SET OF INPUT DATA. IT WAS DESIGNED TO BE 
USED EXCLUSIVELY WITH PROGRAM GECAP. THE TWO MAINLINES ARE 
TIED TOGETHER BY THE USE OF SYSTEM ROUTINE LINK. 


INTEGER P.R 

DIMENSION A ( 12.13) .SUM (20) .WORK (20) . I HLD ( 20 ) .V ( 10) 
DIMENSION Fl<49) *F2(49) 

COMMON TITLXI 7) . T I TLY ( 7 ) *TITLE(7) 

COMMON X ( 101 ) .Y ( 101 ) 

COMMON KI . XMAX . YMAX . X I NC . XST R T 
COMMON YSTRT.XSCLE .YSCLE 
COMMON NCASE . I CASE 
D AT AA / 156*0 . / 


BSTFOO 10 
BSTFOO 20 
BSTF0030 
BSTF 0040 
BSTF0G50 
BSTF0060 
6STF0070 
BSTFU08U 
8STF009U 
BSTF 01 OO 
BSTF01 10 
BSTF 0 1 2U 
BSTF01 30 
BSTF014U 
BSTF015U 
BSTF0160 
BSTF0170 


DATAF 1/161. .18.5 .10.1 .7. 7 1 « 6 . 6 1 * 5 . 99 . 5 . 59 • 5 . 32 » 5 . 12 .4. 96 .4 . 84 .4 • 75BSTF0 180 


1*4.67* 4.60. 4. 54. 4. 49. 4, 45. 4. 41. 4. 38. 4. 35. 4. 32. 4, 30. 4. 28. 4. 26. 4. 24. BSTF0190 
24.23*4.21*4.20.4.18.4.17.4.16*4. 1 5 . 4 . 14 * 4. 1 3 . 4 . 1 2 . 4 . 1 1 , 4. 105 .4. 10 . BSTF0200 
34. 09 *4*08,4.075. 4. 07, 4. 065. 4. 06,4. 05 5,4. 05. 4. 04 5.4. 04. 4. 035/ BSTF0210 

DAT AF2/200. *19.0.9.55 .6.94,5.79.5.14 ,4 . 74 ,4 . 46 . 4 . 26 , h . 1C . 3 .93 . 3 . 89BSTF0220 
1 *3.8 1 *3. 74 » 3.68. 3.63. 3. 5 9. 3. 55. 3. 52. 3.49.3. 47, 3.44, 3.42. 3.40, 3. 39. BSTFC230 


23.37*3.35,3.34*3.3 3,3.32*3.31.3.29,3.28 5 .3.2 8.3.2 7,3.26,3.25 ,3. 24, BSTF0240 
33.235*3.23.3.225,3.22,3.215.3.21.3.205.3.20,3.195.3.19,3.185/ BSTF0250 


C 

C*»* RE-ESTABLISH THE ORIGIN ( SEE MAINLINE GECAP ) 

c 

CALL SCALE(XSCLE.YSCLE.XSTRT.YSTRT) 

KI * KI - 1 
ERROR • 0. 

SAVE 1 - 0. 

SAVE2 *9999999 . 

SAVE3 -9999999. 

NSUB1 ■ 0 
NSUB2 - 0 
NX — 1 
NCOUT-O 
NCT-0 
P - 3 
R * 2 

10 CONTINUE 
C 

C**« BEGIN CURVE-FIT I TER I T I ON 
C 

20 NX - NX + 1 
NY » NX ♦ 1 
C 

C*** LIMIT THE DEGREE OF POLYNOMIAL TO ONE LESS THAN NUMBER OF POINTS 
C**# ON INPUT 


BSTF 0260 
BSTF0270 
BSTF0280 
BSTF0290 
BSTF 0300 
BSTF03 1 0 
BSTF 0320 
BSTF0330 
BSTF0340 
BSTF035 0 
BSTF0360 
BSTF0370 
BSTF 0380 
BSTF0390 
BSTF 0400 
BSTF0410 
BSTF04 2 0 
BSTF 0430 
BSTF 0440 
BSTF045 0 
BSTF0460 
BSTF0470 
BSTF 0480 
BSTF049U 
BSTF0500 


C BSTF0510 

IFINX-KI 130*200,200 BSTF0520 

30 IFINX-9140. 40*200 BSTF0530 

40 CALL CURF(X.Y*KI .NX . NY . A .SUM . V . WORK , I HLD ,E ) BSTF0540 

IF(E)70*70,50 BSTF0550 

50 WR I TE ( P * 60 ) BSTF0560 
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60 FORMAT < / / t ' ***ERROR*** ERROR IN ROUTINE CURF. EXECUTION CONT I NU I BSTF05 70 


1NG') BSTF0580 

70 CONTINUE BSTF0590 

ERROR « 0.0 BSTF0600 

SIGSQ » 0.0 BSTF0610 

C BSTF06 20 

C*** CALCULATE ERRORS FOR THIS POLYNOMIAL AND SUM THEM UP BSTF0630 

C BSTF06 40 




DO 80 KIC « 1.ICI 


BSTF0650 



YCALC « FIA.NX.XOCK) ) 


BSTF066U 



SIGSQ = ABS! Y(KK)-YCALC)**2. 


BSTF0670 


80 

ERROR » ERROR + SIGSQ 


BSTF0680 


90 

CONTINUE 


BSTF0690 


100 

ERROR-ERROR/ (KI-INX+ 1)) 


BSTFU7Q0 



NSUB1 « NX 


BSTF07 10 



NSUB2 ■ NY 


BSTF0720 



SAVE1-SAVE2 


BSTF0730 



SAVE2 « SAVE3 


BSTF0740 



SAVE3 « ERROR 


BSTF0750 


110 

' TEST 1« ABS (SAVE2-ERROR ) /ERROR 


BSTF0760 



TEST 2* ABSI SAVE 1-ERROR ) /ERROR 


BSTF0770 



I-KI-INX+1 ) 


BSTF07 80 



I F ( 1-49 ) 130.130. 120 


BSTF0790 


120 

I«49 


BSTF0800 


130 

I F ( TEST1-F 1 1 I ) 1150.140.140 


BSTF08 10 


140 

NCOUT-O 


BSTF0820 



GO TO 160 


BSTF08 30 


150 

NCOUT * 1 


BSTF0840 


160 

IF ( TEST2-F 2 ( 1 1)180.170.170 


BSTF0850 


170 

NCT-NCOUT 


BSTF0860 



GO TO 20 


BSTF0870 


180 

IFINCT-l 1170.190.190 


BSTF0880 


190 

NSUB1-NX 


BSTF08 90 



NSUB2 ■ NX+ 1 


BSTF0900 



GO TO 220 


BSTF0910 

c 




BSTF09 20 

c#** 

COMPUTE THE BEST-FIT CO-EFFICIENTS 


BSTF09 30 

c 




BSTF0940 


200 

WRITE (P.210) 


BSTFU950 


210 

FORMAT! IX.' ***USERS MESSAGE*** THE 

POL YNOM I NAL CALCULATED 

MAY OR BSTF096Q 


1MAY NOT BE THE BEST FITTED CURVE') 


BSTF0970 


220 

CALL CURFIX.Y.ICI .NSUB1 .NSUB2 . A . SUM . 

V . WORK . I HLu » E ) 

BSTF0980 



WRITE! P.230) 


BSTF0990 


230 

FORMAT! 1H1 ) 


BSTF 1000 



WRITE!P.240)NSUB1 


BSTF1010 


240 

FORMAT!/.' *CO-EFFICIENTS FOR POLYNOMIAL OF DEGREE ' . I 3 . ' * ' 

) BSTF 1020 



DO 260 K«1 .NSU92 


BSTF1030 



J « NSUB2 + 1-IC 


BSTF 1040 



WRITE(P.250)A( J.l ) 


BSTF 1050 


250 

FORMAT ( F20. 12 ) 


BSTF 1060 


260 

CONTINUE 


BSTF 1070 



WRITEIP.270) 


BSTF1080 


270 

FORMAT!/////) 


BSTF 1090 



WRITE(P.280)NSUB1 


BSTF 1100 


260 

FORMAT (• *CALCULAT ED VALUES AND SUM 

OF ERRORS SQUARED FOR' 

.12.' DEBsTFll 10 


1GREE POLYNOMIAL*' ) 

SIGSQ ■ 0.0 

ERROR ■ 0.0 

DO 300 I -1 • K I 

YCALC « FIA.NSUB1 .X! I ) ) 

SIGSQ ■ ABS1YII >-YCALC>**2. 


BSTF1120 

BSTFU30 

BSTF1140 

BSTF1150 

BSTF1160 

BSTFI170 


x. 


52 



n nn 


WRITE IP *290) I .YCALC.SIGSO 


BSTF1180 


290 

300 

3X0 


C 

c*** 

c 


320 

330 
3 AO 


350 

360 

370 


»*# 

380 


FORMAT!' FOR POINT'. 13.' YCALC = ' *F20 . 6 . 1 OX . ' D I FF**2 = '.F20.12) BSTF1X90 


ERROR » ERROR + SIGSQ 
WRITEIP.310IERROR 

FORMAT (//*• SUM OF SQUARES OF ERRORS 
WRITEIP.230) 

plot the calculated function on the 

XBEGN-X ( 1 ) 

STEP ■ XBEGN 
Y I - F( A.NSUB1.STEP) 

CALL EPLOTI l.STEP.YI ) 

STEP * . 1*X I NC 
XLAST *X I K I ) 

X I “XBEGN 
XI “XI+STEP 

I F ( XLAST-X I 13 30. 3 3 0.3 AO 
XI =XLAST 

YI* F I A » NSU81 »X I ) 

CALL EPLOTI2.XI *YI ) 

IF (XLAST-X I 13 50.350.320 
CONTINUE 

IF INCASE- 1 CASE) 380.380.370 
CALL EPLOT(l.XSTRT.YSTRT) 

CALL FINSH 
CALL LI NX I GECAP ) 

SET UP THE PEN FOR THE NEXT PLOT 

CALL EPLOTI 1. XI . 1 1 ) 

CALL EPLOT I l.XMAX.YSTRT ) 

CALL SCALE! 1.0.1. 0.0. 0.0.0) 

CALL EPLOT I 1. A. 0.0.0 ) 

CALL EXIT 
END 


BSTF1200 
BSTF12 10 

FOR CURVE-FIT ='»F20.12) BSTF1220 

BSTF1230 

BSTF12AO 

SUPPLIED GRID I SEE MAINLINE BSTF1250 

8STF1260 
BSTF1270 
BSTF1280 
BSTF1290 
BSTF1300 
BSTF 13 10 
BSTF1320 
BSTF 1330 
BSTF13A0 
BSTF 1 3 50 
BSTF 1360 
BSTF 1 370 
BSTF 1380 
BSTF 1390 
BSTF 1 AOO 
BSTF 1 A 10 
BSTF1A20 
BSTF1A30 
BSTF1AA0 
BSTF1A50 
BSTF 1A60 
BSTF1A70 
BSTF 1A80 
BSTF1A90 
BSTF 1500 
BSTF 15 10 
BSTF 1 520 
BSTF 1530 


// DUP 

•DELETE 

•STORE 


BESFT 
WS UA BESFT 
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// FOR 

♦EXTENDED PRECISION 
♦ONE WORD INTEGERS 

FUNCTION FIA.N.X) FUOOU 

F00 10 

♦ ♦« GECAP REVISION Ct 09/05/72 F0U2u 

FUO 30 
FOOAu 

♦ ♦* FUNCTION F CALCULATES A Y VALUE FRGi-i A GIVEN X VALUE USING THE LEAST FOO&U 

♦♦♦ SOU ARES CO-EFFICIENTS F0u6u 

FU070 

DIMENSION A < 1 2 * 1 3 ) FOOSu 

ID F = A ( N+l * 1 ) FOOVU 

L * N FOiDU 

I F ( N ) 20 * AO » 20 F0110 

20 DO 30 I = l.N F012J 

F = A ( L * 1 ) *X*» I + F F01A0 

30 L - L - 1 F 0 1 5 0 

AO RETURN F0160 

END F0170 

// DUP 

♦DELETE F 

♦STORE WS UA F 
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// FOR 

♦ONE WORD INTEGERS 
♦EXTENDED PRECISION 

SUBROUTINE CURF ( X . Y .N »M .MX . A . SUM . V .WORK . I HLD . E ) 

C 

C*** GECAP REVISION C. 09/05/72 

CATEGORY 

MATHEMATICAL 
PURPOSE 

FIT TABULAR DATA TO POLYNOMIAL OF TYPE Y=A0+A1*X+A2*X**2+. • . 
DESCRIPTION 

THIS SUBROUTINE WILL FIT N DATA POINTS IN X AND Y TO A 
POLYNOMIAL OF THE TYPE Y = AO+Al*X-t-A2»X**2+ • . .+AM*X**M . 

WHERE N IS GREATER THAN OR EQUAL TO M+l. THE COEFFICIENTS 
AO. Alt... .AM ARE DETERMINED BY THE METHOD OF LEAST SQUARES. 
THE ROUTINE MAY BE RUN IN EITHER SINGLE OR DOUBLE PRECISION 
ARITHMETIC. THE DATA POINTS NEED NOT BE EQUALLY SPACED. 
INPUT 

1 X FIRST LOCAT+ON OF AN ARRAY CONTAINING N 

INDEPENDENT VARIABLE DATA POINTS. 

2 Y FIRST LOCATION OF AN ARRAY CONTAINING N 

DEPENDENT VARIABLE DATA POINTS. 

3 N NUMBER OF DATA POINTS TO BE FITTED. 

A M DEGREE OF POLYNOMIAL TO WHICH DATA IS TO BE 

FITTED. 

5 MX VARIABLE DIMENSION INTEGER. MUST BE EQUAL TO 

MAXIMUM M USED PLUS I. 

6 A A TWO DIMENSIONAL ARRAY USED BY CURFIT TO SOLVE 

THE LEAST SQUARES NORMAL EQUATIONS. A MUST BE 
DIMENSIONED TO (MX.MX+l) IN THE CALLING PROGRAM. 

7 SUM A TEMPORARY STORAGE AREA. MUST BE DIMENSIONED 

TO (2*MX) IN CALLING PROGRAM. 

8 V A TEMPORARY STORAGE AREA. MUST BE DIMENSIONED 

TO (MX) IN CALLING PROGRAM. 

9 WORK SEE 8 ABOVE. 

10 I HLD SEE 8 ABOVE. 

OUTPUT 

1 A A(l.l) CONTAINS THE COEFFICIENT AM. A ( 2 • 1 > 

CONTAINS AM-1 .... .AiM+l.l > CONTAINS AO OF THE 
EQUATION Y»A0+A1*X+A2*X**2+. . .+AM*X**M. 

2 E ERROR CHECK. IF E'0.,O.K. IF E=l. AN ERROR 

HAS OCCURED IN THE CALCULATION OF THE 
COEFFICIENTS. 

REMARKS 

THIS ROUTINE CONTAINS A METHOD OF SCALING WHICH PREVENTS 
OVERFLOW IN CASE THE DATA POINTS ARE LARGE. AND PREVENTS 
LOSS OF ACCURACY IN CASE THE INDEPENDENT DATA POINTS 
AGREE TO SEVERAL SIGNIFICANT FIGURES. 

REFERENCES 

INTERNAL TECHNICAL NOTE NO. TN-65-01 
LEAST SQUARES POLYNOMIAL CURVE FIT NO. 2 
BILLY G. GIBBS 

SCIENTIFIC DIGITAL PROGRAMMING BRANCH 
ARMY COMPUTATION CENTER 
REDSTONE ARSENAL. ALABAMA 
CONTACT 

CHIEF. THEORETICAL PROBLEMS SECTION 
SCIENTIFIC DIGITAL PROGRAMMING BRANCH 
ARMY COMPUTATION CENTER 


CURF0000 
CURF0010 
CURF0020 
CURF 0030 
CURF0Q40 
CURF0050 
CURF 0060 
CURF 0070 
CURF00B0 
CURF0090 
CURF0100 
CURF 01 10 
CURF0120 
CURF01 30 
CURF 01 AO 
CURF 0150 
CURF0160 
CURF0170 
CURF0180 
CURF0190 
CURF0200 
CURF0210 
CURF0220 
CURF02 30 
CURF0240 
CURF 02 50 
CURF0260 
CURF0270 
CURF02 80 
CURF0290 
CURF0300 
CURF0310 
CURF0320 
CURF0330 
CURF0340 
CURF0350 
CURF0360 
CURF03 70 
CURF0380 
CURF0390 
CURF0400 
CURF0410 
CURF0420 
CURF0430 
CURF0440 
CURF0450 
CURF 0460 
CURF0470 
CURF0480 
CURF0490 
CURF0500 
CURF05 10 
CURF0520 
CURF 05 30 
CURF 0540 
CURF05 50 
CURF0560 
CURF05 70 
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SUBROUTINE CURE ( X » Y ,N »M .MX » A « SUM » V * WORK. » I HLD * E ) 

CURF0580 


DIMENSION X ( 1 > . Y ( 1 1 »SUM< 1),V(1).A<12,1> .WORK! 1 ) » IHLD(l) 

CURF0590 


E = 0.0 

CURF0600 


LS * 2*M + 1 

CURF06 10 


LB * M + 2 

CURF0620 


LV = M + 1 

CURF 06 30 


XH = ABS < X ( 1 ) ) 

CURF0640 


DO 20 I =2 »N 

CURF 0650 


IFIXH- ABSIXI I) ) 110.20.20 

CURF0660 

10 

XH* ABS ( X ( I ) ) 

CURF 06 70 

20 

CONTINUE 

CURF0680 


DO 30 1*1. N 

CURF06 90 

30 

X< I ) - X < I ) /XH 

CURF0700 


fhh*x ( l ) 

CURF07 10 


DO 50 I * 2 . N 

CURF07 20 


IFIFHH-XI I ) 150.50.40 

CURF0730 

40 

FHH*X ( I 1 

CURF0740 

50 

CONTINUE 

CURF0750 


DO 60 1*1. N 

CURF0760 

60 

X ( I 1 »X ( I l-FHH 

CURF0770 


DO 70 J= 2. LS 

CURF07 80 

70 

SUM ( J ) * 0 . 

CURF0790 


SUM! 11* N 

CURF 0800 


DO 80 J* 1 . LV 

CURF 08 10 

00 

V< J 1*0. 

CURF08 20 


DO 100 1*1. N 

CURF 08 30 


P*l. 

CURF 0840 


V(ll* VI1I+ Y(I) 

CURF08 50 


DO 90 J* 2.LV 

CURF 0860 


P * XI I ) * P 

CURF0870 


SUM ( J ) * SUM ( J ) + P 

CURF 0880 

90 

V(J) * V< J) + Y ( I 1 * P 

CURF 0890 


DO 100 J- LB.LS 

CURF 0900 


P • X< I 1 * P 

CURF09 10 

100 

SUM ( J 1 * SUM ( J ) + P 

CURF 09 20 


KK * LV 

CURF09 30 


DO 120 I * 1 . LV 

CURF0940 


L * I 

CURF 0950 


DO 110 K* 1 » LV 

CURF0960 


A < K .KK ) - SUM ( L 1 

CURF0970 

110 

L * L+ 1 

CURF0980 

120 

KK * KK - 1 

CURF0990 


DO 130 1*1. LV 

CURF1000 

130 

A( I .LSI - V ( I ) 

CURF 1010 


DO 140 1*1, LV 

CURF1020 

140 

I HLD ( I ) *1 

CURF 1030 


JJ*LB 

CURF 1040 


DO 280 1*1. LV 

CURF 10 60 


KK-LV-I 

CURF 1060 


IF(KK)240. 240.150 

CURF1070 

150 

LL*KK+1 

CURF1080 


I JJ*1 

CURF 1090 


L-I 

CURF 1100 


WORK ( 1 ) *A ( 1,1 1 

CURF 1110 


DO 170 I I-l.LL 

CURF 1 1 20 


DO 170 J-l.LL 

CURF 1 1 30 


IF (ABS (WORK! 1 1 l-ABSI A( I I .J ) 1) 160 . 170 , 170 

CURF1140 

1 60 

WORK ( 1 ) * A ( I I »J) 

CURF 1150 


L-J+I-l 

CURF1160 


I JJaJ 

CURF11 70 

170 

CONTINUE 

CURF11 80 
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IF I I JJ-1 1200.200*180 
180 DO 190 I I * 1 * L V 
Z = A( 11,1) 

At II .1 ) =A< I I . I JJ) 

190 At I I , I JJ 1 = Z 
I Y* IHLD ( I 1 
IHLDI I ) "IHLDIL J 
IHLD (LI * I Y 
200 DO 230 L=1.KK 

IFtABSt All ,1 ) l-ABS (AtL + 1 .1) 1)210.230.230 
210 DO 220 J"1,JJ 
Z»A( 1. J) 

At 1 » J ) - A ( L+ 1 , J ) 

220 A<L+1.J)-Z 
230 CONTINUE 
2 AO JJ-JJ-1 

IF (At 1 .1 1 1 250.A30.250 
250 DO 260 J-l.JJ 
280 WORK ( J ) =A ( l.J + 1) /At 1.1) 

KK-JJ+1 
DO 270 K»1,M 
DO 270 J-2.KK 

2 70 A(K.J-l)«A(K+l.J)-A<K+l.l)»WORKl J— 1 > 

DO 280 J-l.JJ 
280 A(LV.J)=WORK< Jl 
LB-LV-1 
DO 310 I-l.LB 
L-I+l 

DO 310 J»L»LV 

I F ( IHLD ( I ) -IHLDI J) ) 310.310,290 
290 I Y* IHLD I I ) 

IHLDI I ) “IHLDI J) 

IHLDI J) =IY 
DO 300 K-l.l 
Z«A( I .to 
All ,K ) -At J.K ) 

300 A ( J , K ) -Z 
310 CONTINUE 

DO 320 I «1 ,N 
321 X ( I ) "X ( I H-FHH 
NN-LV 

DO 380 I "1 « LV 
SUM ( 1 ) "0. 0 
NN-NN-1 
KK-NN+1 

IF(KK-M)330»330»380 
330 L-M-XK+1 

DO 370 K-KK.M 
I I "X-NN 
SIGN-1. 

FAC-1. 

I F ( I 1-0) 3A0.360.3A0 
3 AO SIGN" ( -1 ) ** 1 1 
IF-K 

DO 350 J"1,II 

Z"J 

WORK 1 1 ) * I F 
FAC"FAC*WOR<( 1)/Z 
350 IF" IF-1 

360 SUM ( 1 ) "SUM ( 1 ) +FAC*SI GN*FHH**I I *A ( L , 1 ) 

370 L-L-l 


CURF 1190 
CURF 1200 
CURF1210 
CURF 12 20 
CURF 12 30 
CURF 12 AO 
CURF12 50 
CURF1260 
CURF1270 
CURF1280 
CURF 1290 
CURF1300 
CURF 1 310 
CURF 13 20 
CURF1330 
CURF13A0 
CURF1350 
CURF 1360 
CURF1370 
CURF 1380 
CURF 1390 
CURF1A00 
CURF1A10 
CURF1A20 
CURF1A30 
CURF 1 A AO 
CURF1A50 
CURF 1 A60 
CURF1A70 
CURF1A80 
CURF1A90 
CURF1500 
CURF 1510 
CURF 1520 
CURF1530 
CURF15A0 
CURF1550 
CURF1560 
CURF 1570 
CURF1580 
CURF 1590 
CURF 1600 
CURF 1610 
CURF 16 20 
CURF 16 30 
CURF 16 AO 
CURF 1650 
CURF 1660 
CURF 1670 
CURF 16 80 
CURF 1690 
CURF17 00 
CURF 1710 
CURF 17 20 
CURF 1730 
CURF 17 AO 
CURF1750 
CURF1760 
CURF 1770 
CURF 1780 
CURF1790 
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380 A ( I .2>=A< I t 1 ) +SUM ( 1 ) 

DO 390 1*1 * LV 
390 A( I . 1 ) “A ( I ,2) 

DO AOO I * 1 * N 
A 00 X( I )*X( I )*XH 
K«LV+1 

DO A10 I » 1 * L V 

K*K-1 

NN«I-1 

A10 A ( K • 1 ) ■ A ( K t 1 ) / ( XH**NN ) 
A20 RETURN 
A30 E-l. 

GO TO A20 
END 

// DUP 

•DELETE CURF 

•STORE WS UA CURF 


CURF 1800 
CURF 1810 
CURF 1820 
CURF 1830 
CURF 18 AO 
CURF 1850 
CURF 1860 
CURF 1870 
CURF 1880 ' 
CURF 1890 
CURF 1900 
CURF 1910 
CURF 19 20 
CURF 1930 
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1 


// ASM 


ENT FINSH 

FINSH DC 0 

MDX L 50*0 
MDX *-3 

BSC I FINSH 
END 


// DUP 

•STORE WS UA FINSH 


ASM 1 

ASM 

ASM 3 

ASM A 

ASM t> 

ASM 6 

ASM 7 

ASM fa 
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APPENDIX C 


SYMBOLS AVAILABLE FOR THE IBM 1130 PLOTTER ROUTINES 


Alphabetic Characters 
Numeric Characters 
Special Characters 
Characters 


( 

+ 

$ 

* 

) 

/ 


Space 


A - Z 
0-9 


Punch 

12 - 8-3 

12 - 8-5 

12 - 8-6 

11 - 8-3 

11 - 8-4 

11 - 8-5 

11 

0-1 

0 - 8-3 

8-5 

8-6 

Blank 
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APPROVAL 


COMPUTER USER'S MANUAL FOR A GENERALIZED 
CURVE FIT AND PLOTTING PROGRAM 


By Ronald A. Schlagheck, B. D. Beadle II, 

B. D. Dolerhie, Jr., and J. W. Owen 

The information in this report has been reviewed for security classifi- 
cation. Review of any information concerning Department of Defense or 
Atomic Energy Commission programs has been made by the MSFC Security 
Classification Officer. This report, in its entirety, has been determined to 
be unclassified. 

This document has been reviewed and approved for technical 
accuracy. 




Chief, Mechanical and Crew Systems 
Integration Division 


Uuvi L- 


KARL L. HEIMBURG 

Director, Astronautics Laboratory 
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